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A  REVIEW  OF  REPRESENTATION  FUNCTIONS 


FOR  PROBABILITY  OF  DETECTION 

by 

Wayne  Rivers,  J.  N.  Bucknam, 
E.  N.  Khoury,  and  R.  E.  Blase 


ABSTRACT 


Simple  representation  functions  that  interrelate  the  primary 
signal  detection  variables  for  the  receiver  structure  of  Marcum  and 
Swerling  are  reviewed  in  regard  to  accuri  .  complexity  and  inversion. 
Functions  reviewed  include  those  of  Brooks,  Neuvy,  and  Khoury-Bucknam. 

The  first  two  of  these  are  simple  algorithms  for  computing  minimum  ratio 
of  signal  energy  to  noise  power  density  as  functions  of  number  of  samples 
integrated,  the  target  distribution  case,  and  required  probabilities  of 
detection  and  false  alarm.  The  Khoury-Bucknam  functions  are  analytic 
and  invertable,  and  they  relate  probability  of  detection  and  signal-to- 
noise  ratio  using  parameters  chosen  uniquely  for  each  case,  number  of 
samples,  and  probability  of  false  alarm. 

The  procedures  and  software  that  support  determination  of  the 
coefficients  of  the  Khoury-Bucknam  function  are  documented. 
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1.  INTRODUCTION 


The  process  of  determining  when  a  desired  signal  is  present 
in  a  radar  receiver,  called  detection,  was  analyzed  by  Marcum  and  Swerling 
[1]*  in  terms  of  the  receiver  structure  shown  in  Figure  1.  In  the  con¬ 
text  of  Marcum's  and  Swerling' s  analysis,  the  detection  process  is  de- 

scribable  in  terms  of  functions  of  five  (total)  variables: 

•  A  target  distribution  identifier  and  its  associated  number 

of  degrees  of  freedom  in  the  target  signal  (K) ; 

•  The  number  of  signal  samples  integrated  at  video  (N) ; 

•  The  probability  that  the  video  integrator  output  exceeds 

a  threshold  level,  given  that  no  target  signal  is  present: 

the  "probability  of  false  alarm"  (P_,); 

r  A 

•  The  probability  that  the  video  integrator  output  exceeds 

the  threshold  level,  given  that  the  target  signal  is  present: 
the  "probability  of  detection1'  (PD) ; 

•  The  ratio  of  energy  in  a  sample  of  the  "target”  signal 

to  the  background  noise  power  density,  the  so-called  "signal- 
to-noise  ratio"  (S) . 

In  addition  to  the  stipulated  receiver  structure,  other  parameters  of  the 
problem  were  fixed:  Che  Gaussian  distribution  of  the  background  noise 
with  known  fixed  average  power  (or  energy  density) ;  ideal  and  pure  target 
distribution  cases  belonging  to  the  chi-square  family;  known  time  of 
arrival  of  the  target  signal  (but  not  phase);  ideal  device  linearity, 
or  in  the  case  of  the  envelope  "detector"  ideal  square-law  transfer 
properties;  perfect  matched  filter  on  the  coherent  side;  large  video 
bandwidth  on  the  video  side  compared  to  the  matched  filter  bandwidth; 


Numbers  in  brackets  refer  to  references  in  Section  4. 
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satisfaction  of  the  narrow-band  assumption  for  the  coherent  matched  fil¬ 
ter;  and  perhaps  others. 


For  a  given  set  of  radar  requirements  or  an  existing  radar, 
appropriate  values  will  be  known  for  the  number  of  samples  integrated, 
the  probability  of  false  alarm  and  the  target  distribution.  Therefore, 
the  last  two  variables  above  stand  out  as  principal  candidates  for  in¬ 
dependent  and  dependent  status.  In  modeling  radar  performance,  compu¬ 
tation  of  probability  of  detection  for  values  of  signal-to-noise  ratio 
is  a  routine  and  highly  repeated  activity.  The  inverse  of  that  function, 
the  signal-to-noise  ratio  needed  to  achieve  a  given  probability  of  de¬ 
tection,  is  also  widely  used,  but  in  a  less  repetitious  mode. 

4 

Marcum  and  Swerling  [1]  published  functions  of  PD(S=(R/R^)  ) 

as  plotted  curves,  with  P_.,  N,  and  Case  as  parameters,  and  the  signal- 

r  A 

to-noise  power  ratio  expressed  indirectly  in  terms  of  a  normalized  range 

variable.  Meyer  and  Mayer  [2]  made  computations  and  presented  results 

as  plotted  curves  of  P_(N),  with  ,  n'  =  ln(2)/F_.,  Case  as  parameters. 

u  Db  r  A 

This  latter  set  of  curves  is  eminently  readable  and  suitable  for  manual 
methods  of  model  computation,  but  the  procedure  for  their  computation 
by  machine  is  lengthy.  Fehlner  [3]  implemented  a  series-representation 
of  P_(S,  N,  Y,  (P_,  .Case,  N) ,  Case)*  which  is  in  use  today  as  part  of  the  NRL 
model  package  as  subroutine  MARSWR  [e.g.  4).  Computation  cost  using 
this  routine  can  be  high  (30  to  50%  of  the  total  in  generating  plots  of 
detection  performance)  and  inversion  to  compute  S(PD)  is  possible  only 
by  iteration.  Mitchell  and  Walker  [5]  defined  a  simpler  more  compact 
algorithm  than  MARSWR  in  terms  of  recursive  properties  of  the  functions 
that  are  evaluated  by  the  long  series  in  the  other  methods.  However, 
its  computation  costs  can  be  even  higher  [5]. 


Y  (P  ,  Case,  N)  is  the  threshold  level  in  video  voltage  units, 
b  i'  A 


Over  the  years,  various  approximation  methods  have  been  defined 

to  reduce  the  complexity  and  cost  of  computation,  at  the  expense  of  accuracy. 

North  [6]  assumed  that  signal-to-noise  ratios  of  most  interest  are  high 

so  that  an  asymptotic  form  of  the  modified  Bessel  function  could  be  used, 

resulting  in  a  marked  simplification  of  the  probability  integral  at  a 

cost  in  accuracy  of  only  0.2  dB.  Brooks  [7]  fit  the  North  approximation 

to  produce  a  simple  algorithm  amenable  to  pocket-calculator  implementation 

of  the  function  S(PD>  N,  Pp^,  Case).  This  algorithm  is  examined  further 

in  the  next  section.  Neuvv  [8]  under  Skolnik's  guidance  factored  the 

function  S(Pn,  N,  P_.  ,  Case)  into  the  product  S  =  f.(P_.)  f~(N,  Case) 

D  FA  1  FA  2 

f^Pp,  <'ase)  *  an<*  the  three  factors  with  simple  functions.  This 

method  is  also  examined  in  Section  2. 

Khoury  and  others  at  TSC  recognized  the  need  for  a  representation 

both  efficient  and  reasonably  accurate  for  the  function  r^vS)  which  also 

possessed  ar.  inverse  S(PD).  Khoury  identified  such  a  function  having 

four  parameters:  A,  B,  C,  and  P„,: 

r  A 

PD  -  1  -  (1  -  PpA)/(l  +  (A-S)B)C,  (1) 


and  its  inverse 


(2) 


Bucknam  defined  S  *  10 
Dd 


log^g(S)  and  rewrote  the  functions  as 


1  -  1-P 


FA 


1+10 


0.1(SDB+  AD)\B\C 


(3) 
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1/C 


-  AD 


(4) 


in  which  AD  -  10  log  g(A).  Formula  4  is  linear  in  the  coefficients  AD 
and  1/B,  which  fact  has  significance  for  the  procedure  which  defines  those 
parameters . 
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Khoury  used  sets  of  pairs  of  (Pq»  S)  Input  to  an  iterative  pro¬ 
cedure  for  finding  values  of  A,  B,  and  C  that,  apply  to  a  particular  para¬ 
meter  set  (Case,  N,  P^)  of  interest.  This  procedure  is  the  subject  of 
Sections  3.1  and  3.2  of  this  report.  Buckuam  derived  a  procedure  taking 
advantage  of  the  asymptotic  properties  of  Formula  A  and  its  linear  form, 
and  this  procedure  is  the  subject  of  Section  3.3  of  this  report. 


2 .  COMPARISON  OF  REPRESENTATION  ACCURACIES 


t* 

i 

i 


The  use  of  approximating  representations  is  appealing  because 

of  the  much  reduced  complexity  and  computation  effort.  But  before 

approximations  can  be  accepted,  users  must  appreciate  the  consequences 

in  terms  of  the  accuracy  of  representation.  Here  the  Brooks,  Neuvy, 

and  Khoury-Bucknam  (K-B)  representations  are  compared  with  values  read 

from  the  Mever  and  Mayer  curves  [2].  Values  of  Sno  are  tabulated  for 

Do 

all  four  methods  over  the  following  space: 


PD  -  0.01,  0.1,  0.32,  0.5,  0.8,  0.9,  0.99 


P  =  10“6,  10“8,  10~10 


.  N  =  1,  10,  100 

Case  =0,  1,  2,  3,  4 

and  the  results  are  shown  in  Tables  i-5.  Computation  using  Brook's  and 
Neuvy 's  algorithms  was  straightforward  and  in  accord  with  the  formulas 
in  Appendix  A  of  this  report.  Coefficients  AD,  B,  and  C  were  obtained 
using  Bucknam's  Quick  PDFIT  recipe  applied  to  the  Meyer  and  Mayer  data 
set  as  input,  and  then  Formula  4  was  used  to  compute  the  K-B  entries. 

The  K-B  function  using  Quick  PDFIT  coefficients  departs  from 
the  M&M  values  least  for  Case  0  (non-fluctuating  target)  and  greatest 
!  for  Case  1  (Correlated  Rayleigh  target),  with  Cases  4,  2,  and  3  inter¬ 
im  mediate.  For  Cases  0  and  2  (N  _>  10)  deviations  of  less  chan  0.2  dB  are 

t  typical,  whereas  for  Case  1  0.6  dB  differences  occur  over  0.1  <_  <_  0.9 

and  1  dB  over  0.01  <_  <_  0.99.  The  0.2  dB  error  estimate  is  probably 

j  not  significant,  because  the  internal  consistency  of  the  M6M  curves  is 

r  about  0.2  dB  peak- to-peak. 

! 


L 


j 

| 


J 


Brook's  formulas  work  well  for  all  cases,  with  errors  generally 

less  than  0.5  dB  for  0.1  <_  <_  0.9  and  less  than  1  dB  for  0.01  £  ^  0.99. 

Neuvy’s  algorithm  fails  badly  for  Case  0,  with  differences  ex¬ 
ceeding  3  dB  at  low  values  of  P^.  For  other  cases,  however,  its  perfor¬ 

mance  is  comparable  to  but  a  little  poorer  overall  than  Brooks'. 

For  record  purposes,  values  of  the  coefficients  A,  AD,  B,  and 
C  obtained  using  Quick  PUFIT,  and  which  were  used  in  this  comparison, 
are  tabulated  in  Appendix  B  of  this  report. 
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3.  DETERMINATION  OF  COEFFICIENTS  OF  THE 


KHOURY-BUCKNAM  FORMULAS 


Here  are  discussed  the  three  principal  tools  currently  in  use 
for  finding  values  of  the  coefficients  A,  B,  and  C,  (also  denoted  as 
PI,  F2,  and  P3)  in  Formulas  1-4  of  Section  1. 

3.1  Program  PDFIT 

PDFIT  calculates  a  set  of  values  of  S  corresponding  to  a  set 
of  12  values  of  P^  using  the  recursive  procedure  of  Mitchell  and  Walker  [5] 
in  an  iterated  successive  approximation  mode.  Using  trial  values  of 
the  three  coefficients  and  Formula  (1)  of  Section  1,  the  mean  square 
error  in  Is  computed.  Each  coefficient  in  turn  is  perturbed  by  a 
small  amount  to  compute  partial  derivatives  and  the  matrix  procedure 
outlined  in  Khoury's  memorandum  in  Appendix  C  of  this  report  is  used  to 
compute  new  values  of  the  coefficients.  This  procedure  is  iterated 
until  either  a  time  limit  is  exceeded  or  the  coefficient  values  have 
settled.  The  input  data  set  (S,  P^) ,  the  computed  values  of  using 

Formula  1,  and  the  error  are  printed  to  the  operator  along  with  the  values 
of  the  coefficients.  Note  that  PDFIT  yields  coefficients  that  minimize 
the  mean  square  error  in  PD(S). 

Figure  3-lA  shows  the  dialogue  between  the  operator  and  the 
program  PDFIT  in  a  typical  example,  and  Figure  3-1B  shows  the  formatted 
output  that  resulted.  The  program  PDFIT  along  with  its  subroutines, 
including  implementations  of  the  matrix  coefficient  adjusting  procedure 
and  the  Mitchell-Walker  procedure  are  given  in  Appendix  D  of  this  report. 


OK,  PDFIT 
GO 


PROGRAM  PDFIT  COMPUTES  AN  EMPIRICAL  FIT  TO  THE 
EXACT  PD  VS.  SNR  CURVE  DESIRED.  THE  EXACT  CURVE  IS 
SPECIFIED  BY  SWERLINO  CASE  NUMBER,  PROBABLITY  OF 
FALSE  ALARM,  AND  NUMBER  OF  PULSES  NON-COHERENTLY 
INTEGRATED.  THIS  PROGRAM  WILL  RETURN  THE  THREE 
PARAMETERS  WHICH  GIVE  THE  BEST  FIT  IN  THE  MMSE 
SENSE  <CF.  REF.  TSC-W21-40) 

PDFIT  WAS  DEVELOPED  BY  E.  N.  KHOURY,  AND  R.  F. 
PUSATER I  -  JUNE  79. 


INPUT  DATA  CONSISTS  OF  THE  FOLLOWING: 

A)  OUTPUT  FILE  NAME 

B )  SWERLINO  CASE  NUMBER 

C)  PROSAALITY  OF  FALSE  ALARM 
D>  NUMBER  OF  CASES  TO  BE  RUN 

E)  NUMBER  OF  PULSES  INTEGRATED  FOR  EACH  CASE 


ENTER  OUTPUT  FILE  NAME  -  6 
OUTW3 


ENTER  SWERLINO  CASE  NUMBER  -  0  1  2 


ENTER  PROBABLITY  OF  FALSE  ALARM 
RANGE:  1.  E-3  THROUGH  1.  E-15 
.  IE-9 


ENTER  NUMBER  OF  CASES  TO  BE  RUN 
RANGE:  1  THROUGH  100 


ENTER  NUMBER  OF  INTEGRATED  PULSES  FOR  EACH  CASE 
RANGE:  1  THROUGH  100 


PROCESSING  CASE  NUMBER:  1 


Figure  3-1A  -  Operacor-Program  Dialogue  in  Use  of 
Program  PDFIT 
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CASE  NUMBER  1 


INPUT 

PARAMETERS 

9WERLIN0  CASE  NO. 

PFA 

PULSES  INTEGRATED 

0.  0 

0. 

IE-03 

to 

PARAMETER  ESTIMATES 


PARAMETER  1 


PARAMETER  2  PARAMETER  3 


0. 3353940 


4.  1603013 


2.  3347212 


FtT 


N 

SNR 

PD 

PDP 

1 

0 

1227E 

01 

0.  04918 

0.  03413 

2 

0 

1443E 

01 

0.  09939 

0.  10280 

3 

0 

1710E 

01 

0.  19304 

r  l‘"?76 

4 

0 

1943E 

01 

0.  30229 

r  ’7 

9 

0. 

21S8E 

01 

0.  39763 

C  J . , /3 

6 

0. 

2304E 

01 

0.  49248 

0.  49093 

7 

0 

2323E 

01 

0.  60708 

0. 60814 

9 

0 

2740E 

01 

0.  70633 

0. 70943 

9 

0. 

2973E 

01 

0.  79379 

0. 79873 

10 

0 

3399E 

01 

0.  90682 

0.  90608 

1 1 

0 

37496 

01 

0.  93613 

0.  93240 

12 

0. 

4803E 

01 

0.  99722 

0.  99403 

ERROR 

-0  00494 
-0.  00342 
0.  00028 
O.  00282 
0.  00290 
0. 00133 
-0.  00106 
-O.  00280 
-O.  00296 
O-  00074 
0.  00373 
0.  00318 


MEAN  ERROR  »  «0.  0000 


RMS  ERROR  -  0.  00283 


Figure  P.-13 


Content  of  Output  File  0UTW3, 


Resulting 


from  Program  PDFIT 
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3.2 


Program  QFIT 

OFIT  operates  on  arbitrary  sets  of  (P  ,  s_n)  data  input  via  a 

previously  prepared  file.  Two  options  are  available  for  deriving  the 

coefficients  of  the  fit  equation.  One  applies  the  Quick  PDFIT  procedure 

to  be  described  in  Section  3.3  and  prints  a  summary  of  the  input  data, 

the  errors,  and  the  coefficients.  This  option  yields  coefficients  that 

minimize  errors  in  S__  computed  using  Formula  (A)  of  Section  1.  The 

Do 

other  option  uses  the  procedure  of  Program  PDFIT  operating  on  the  arbitrary 
input  data  and  produces  the  same  summary  as  PDFIT.  This  option  yields 
coefficients  that  minimize  errors  in  using  Formula  (1)  or  (3)  of 
Section  1.  For  both  options,  the  input  data  and  computed  values  from 
the  fit  formula  can  be  plotted. 

Figure  3-2A  shows  a  sample  dialogue  in  which  both  options 
and  the  plot  are  executed.  Figure  3-2B  shows  the  input  file  for  the 
example  and  Figure  3-2C  the  Quick  PDFIT  output  summary.  Figure  3-2L' 
shows  the  output  summary  for  the  PDFIT  option,  and  the  plot  of  input 
data  with  computed  values  using  both  sets  of  coefficients  is  shown  in 
Figure  3-2E. 

Program  QFIT  Including  all  of  its  subroutines  is  included 
as  Appendix  E  of  this  report. 
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OK,  R  «<3FIT  1 

GO 

ENTER  INPUT  FILE  NAME, FORMAT  6A2 
TEST.  WI8 

TO  USE  BOOS  VERSION  ENTER  l.ELSE  0 
1 

ENTER  OUTPUT  FILE  NAME,  FORMAT  6A2 
OUT. WIB 1 

IF  A  PLOT  OF  RESULTS  IS  DESIRED  ENTER  1,  ELSE  0 
1 

ENTER  XDIM(IN), YDIM(IN) 

7,  3 

ENTER  X  AXIS  LABEL,  FORMAT  32A1 

S/N 

ENTER  YAXIS  LABEL,  FORMAT  32A1 
PD 

TO  DO  COMPLETE  PDF IT  EN1ER  1,  ELSE  0 
1 

PDFIT:  TERMINATION  3.  MAX.  ITERATIONS 


ENTER  OUTPUT  FILE  NAME  -  12 
OUT. WIB2 

ENTER:  CASE  PFA  NP  NB 
0,  .  IE-3,  10,  4 

IF  A  PLOT  OF  RESULTS  IS  DESIRED  ENTER  1.  ELSE  0 
1 

ENTER  XDIM(IN),  YDIM<IN) 

7,  3 

ENTER  X  AXIS  LABEL,  FORMAT  32A1 
ENTER  YAXIS  LABEL,  FORMAT  32A1 


OK,  COMO  -E 


Figure  3-2A  -  Operator-Program  Dialogue  in  Use  of  Program  QFIT 


6 

?.  8,  .  06 
10.  199, .  128 

10.  399,  .  276 
11,.  604 

11.  398,  •  909 
11.  801,  .  981 


Figure  3-2B  -  Content  of  Input  File  Test.WlE,  Manuailv-Fed  Data 


RESULTS  OF  QFIT 


SNOB 

SN 

P 

PFIT 

ERR 

9'. 

8000 

9.  8499 

0.  0600 

0.  0838 

0. 

0068 

10. 

1990 

10. 4689 

0.  12B0 

0.  1388 

-0. 

0078 

10. 

8990 

11.  4789 

0.  2760 

0.  3132 

-0. 

0392 

11. 

OOOO 

12.  8892 

0.  6040 

0.  6092 

-0. 

0082 

11. 

3980 

13.  7978 

0.  9090 

0.  8777 

0. 

0313 

11. 

8010 

IS.  1391 

0.  9810 

0.  9833 

-0. 

0023 

FIT  PARAMETERS 


P-l-< l+{A*SN)**a>**(-C) 


0. 71290E-01 


B 

0.  10780E  02 


0.  34S40E  01 


Figure  3-2C  -  Content  of  Output  File  OUT.WIBl,  Results  of  Fit 
to  Manually-Fed  Data 


H1EBULL  -  .  33 

INPUT  PARAMETERS 

SV)eWLIW0  CA8E  W.  PFA  °UL8E8  INTEORATED 

0-  0  0.  IE-09  10 

PARAMETER  ESTIMATES 


PARAMETER  1 

PARAMETER  3 

PARAMETER  3 

0.  0493103 

10.  9946313 

143.  7387998 

FIT 

N 

SNR 

PD 

PDF 

ERROR 

** 

•• 

1 

0.  9990E  01 

0.  06000 

0.  04773 

0.  01337 

3 

0.  1047E  03 

0.  12000 

0.  13138 

0.  00663 

3 

0.  I148E  03 

0.  27600 

0.  29043 

— 0  01443 

4 

0.  1399E  03 

0.  60400 

0  99776 

0  00634 

9 

0.  1390E  03 

0.  90900 

0.  90896 

0.  00044 

4 

0.  1314E  03 

0.  98100 

0  99817 

-0.  01717 

MEAN  ERROR  »  -0  0010  RMS  ERROR  ■  0.  01103 

Figure  3-2D  -  Content  of  Output  File  0UT.WIB2,  Results  of  Fit 
to  P^,  SNDB  Data  from  Mitchell  &  Walker's  Routine 
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3. 3  The  Quick  PDFIT  Procedure 

Quick  PDFIT  cakes  advantage  of  two  properties  of  Formula  4: 

SDB  '  f  Uho  |  (pff)  *  lj  -  “ 

One  is  its  asymptotic  behavior  near  PQ  ■  0  and  P  ■  1;  the  other  is  its 

linear  form  in  coefficients  AD  and  1/B.  The  first  property  is  used  to 

obcain  an  accurate  estimate  of  the  coefficient  C,  by  operating  on  pairs 

of  (P,.,  Sn„)  values  only  at  the  extremes.  The  second  property  is  then 

used  with  a  simple  least-squares  procedure  operating  on  all  of  the 

(P  ,  S_n)  data  to  estimate  coefficients  AD  and  B  that  minimize  errors 
D  DB 

in  S__  computed  by  Formula  (4) . 

Dd 

Figure  3-3  lists  the  steps  in  the  Quick  PDFIT  algorithm,  and 

Figure  3-4  lists  the  FORTRAN  subroutine  that  implements  it.  Bucknam's 

memorandum  deriving  the  procedure  is  included  in  Appendix  F.  Figures 

3-3  and  3-4  and  Appendix  F  assume  that  PFA  is  either  zero  or  insignificant 

For  finite  P_A,  the  same  procedure  can  be  used  by  fitting  a  distorted 
f  A 

set  of  PD  values,  given  by 


Use  of  the  Quick  PDFIT  procedure  on  the  data  pairs  (P.T,  S  ),  with  P„. 

u  Do  r  A 

assumed  zero,  yields  values  of  AD,  B,  and  C  for  use  in  the  more  general 
expression  of  equation  (4). 

No  claim  is  made  that  the  value  cf  coefficient  C  obtained  by 

the  Quick  PDFIT  procedure  is  optimum  in  any  sense,  although  given  the 

initial  estimate  of  C,  the  other  two  parameters  are  optimum  for  minimum 

errors  in  S^.,.  However,  manual  adjustment  of  values  of  C  about  the  esti- 
Do 

mated  one,  followed  by  recomputation  of  AD  and  B ,  hac  not  been  able  to 
show  any  significant  improvement  in  error  over  that  which  results  from 
the  Quick  procedure.  It  has  been  concluded  that  errors  in  fitting  input 
data,  such  as  seen  in  Section  2  of  this  report,  result  from  intrinsic 
shape  limitations  of  Formulas  (1)  -  (4)  and  not  from  a  poor  estimate  of 
parameter  C. 


(4) 
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Figure  3-3  -  The  Quick  PDFIT  Algorithm 


1.  Define  function  G(X) 


2.  Input  pairs  of  the  variables  SDBi,  i  -  1.  ...  n,  in  order  of 
increasing  values. 

3l  Yi  =  log10  log10(1  *  Vj*  1  “  !•  2 


Y2  ~  Y1 
SDB2  ”  SDB1 


zi  m  ~  loS10(1  ~  i  ■  n-1,  n 


Z  -  Z  , 

M  .  _n _ nil _ 

SDBn  “  SDBn-l 
M 

C  = 


4.  Q.  =  1  -  (l  -  Pi)1/C,  i  =  !,  ...  n 

G(Q  )  -  G(Q  ) 

5,  M  - - 2 _ n-1 

^  I°810(Qn)  -  1°8lO<(J„-l) 

c(Q.)  -  G(Q,) 

Mg  =  - - - - ± - - - 

logio(c  lnV  -  logl0(^  ln(V 

C  =  C  Mj^/Mg  (this  step  adjusts  the  first  estimate  of  C  to  an  improved  value). 
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Figure  3-3  (continued) 


6-  Wi  -  GCP^,  i  -  l,  ...  n 

7.  Uw  ■  EW^,  i  -  1,  . . .  n 


US  “  ZSDBi’  1  "  1»  "•  n 
USS“  n 


USW"  "SDBi  Wi ’  1  “  1 >  • • •  n 


8.  D  -  n  •  Uss  -  (USV 


H1  =  (U3S  “w  -  US  USW)/D 


H2  =  (N  '  USW  -  US  V/D 


9.  AD  = 


B  -  10  H, 


nno  non  non  nno  non  non  nnnnnnnnnnn 


Figure  3-4  Subroutine  QPDFIT 


SUBROUTINE  QPDFIT  <PD,SDB,N,  P> 

COMPUTES  APPROXIMATE  VALUES  FOR  PARAMETERS 
PJ.P2,  P3  IN  THE  EXPRESSION: 

PD«<  l-<  1  +  <PI*S)*-*P2)**-P3)*100 

WHERE:  PD=DETECTIQN  PROBABILITY  (IN  PERCENT) 
S=S/N=AL0G10(SDB/10.  0) 

SEE  TSC-  W7-  73 

DIMENSION  Y(2) ,  ZC2). W(200>,  PD< 1).  P<3), SDB < 1 ) 
DIMENSION  PDD < 200 ) 

REAL  MY,  M2 
DO  100  1*1,  N 

PDD ( I )=PD< I )/100. 0 
100  CONTINUE 

COMPUTE  SLOPE  3ETWEEN  FIRST  TWO  POINTS 
DO  150  1  =  1,2 

Y ( I )»AL0C10(-AL0010( 1.  -PDD ( I ) ) ) 

150  CONTINUE 

MY*  <  Y ( 2 ) — Y ( 1 ) ) / ( SDB  <  2 ) -SDB ( 1 ) ) 

COMPUTE  SLOPE  BETWEEN  LAST  TWO  POINTS 

NM1-N-1 

IZ-0 

DO  200  I-NM1, N 
IZ-IZ+1 

Z( IZ)»-AL0G10< 1. -PDD<  I  >  ) 

200  CONTINUE 

MZ* ( Z ( 2 ) -Z ( 1 ) ) / ( SDB ( N  > -SDB  <  NM1 ) ) 

COMPUTE  INITIAL  ESTIMATE  OF  K3 
P ( 3 ) *MZ/MY 

DISTORT  THE  PD  VALUES 
DO  225  1-1,  N 

PDD ( I > *1 .  -( < 1.  -PDDC I  >>#■*<  1.  /P (3  >  >  > 

225  CONTINUE 

COMPUTE  LARGE  PD  SLOPE  RATIO 

SRL1-AL0G10C  (  <  l.  -PDD<N>  >*•»<-;.  /P(3)  )  )-l.  ) 
SRL2*AL0010( (Cl. -PDD(NM1 >>**(- 1 .  /P  < 3 >) >-l .  ) 

SRL3*( -1.  /Pi  3) ) *ALOO 1 0  < 1.  -PDD(N) ) 

SRL4»(  1.  /P  ( 3  >  >  *ALQG  1 0 <  1.  -PDD ( NM1 )  ) 

SRL*(SRL1-SRL2) /(SRL3+SRL4) 

COMPUTE  THE  SMALL  PD  SLOPE  RATIO 

SRS1»AL0010( ( ( 1.  -PDD ( 2 ) >**(-l.  /P(3) ) >-l.  ) 
SRS2=*AL0G10(  (  ( 1,  -PDDC  1  )  )**(-l.  /P<3>  )  )-l.  ) 
SRS3-AL0Q10<  (-ALOOC  10.  0)  /P(3>  )*<AL0G10<  1.  -PDD<2)  )  )  ) 
SRS4-AL0010< ( — ALOG ( 10.  0) /P(3) >*<AL0G10< 1.  -PDD< 1 ) ) ) ) 
SRS* ( SRS1 -SRS2 ) / ( SRS3-SRS4 ) 

P (3)-P  <3)*(SRL/SRS> 
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c 

c 


c 

c 

c 


290 

C 


C 

C 

C 


C 


initialize  counters  for  linear  regression 

SUM»0 

SUMX-O.  0 

SUMY-O.  0 

SUMX2-0.  0 

SUMXY-O.  0 

COMPUTE  U<1>  AND  PARAMETERS  FOR  A  STRAIGHT  LINE  FIT 
DO  250  1-i.N 

SU^SUmS®1?^1*"  01*PD<  I  ))**<-!.  /P(3)))-i.  > 
SUMX=SUMX+SDB< I ) 

SUMY-SUMY+W<I> 

SUMX2-SUMX2+SDB  < I ) *SDB ( I ) 

CONTI^’<V'SUMXY’SDB<i!*«<^ 

D-SUM*SUMX2-SUMX*SUMX 
B— <SUMX2*SUMY— SUMX*SUMXY ) /D 
A* <  SUMXY*SUM-SUMX*SUMY ) /D 

COMPUTE  K1  AND  K2 


RETURN 

END 
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APPENDIX  A 


Approximation  Algorithms  of  Brooks  and  Neuvy 


Table  A-l  Brooks'  Algorithm  for  Minimum  Required 
Signal-to-Noise  Ratio  [7] 


Reference :  L.  W.  Brooks,  Jr.,  "Detection  Theory  on  a  Pocket  Calculator," 

Proceedings  of  IEEE  Southeastcon  ,  April  197  7;  p  354-357. 


Definitions : 

■  Probability  detection  occurs,  given  that  target  signal 
is  present. 

PpA  -  Probability  that  detection  occurs,  given  that  target 

signal  is  not  present  -*■  probability  of  false  alarm  per 
cell,  per  N-sample  group. 

N  a  Effective  number  of  samples  integrated  at  video 

£  =  Asymptotic  efficiency  of  the  envelope  detector 

K  =  Number  of  degrees  of  freedom  in  the  target  signal  in 
N  samples 

S  =*  Ratio  of  target  signal  energy  to  noise  power  density 
that  would  produce  P^. 

Algorithm: 


Detector 

Square  Law 

1.0 

Linear  Law 

0.915 

Logarithmic  Law 

0.608 

2 .  Target  Distribution  Case _ K_ 

0.  Nonfluctuating  00 

1.  Rayleigh;  Correlated  over  N  1 

2.  Rayleigh;  Fluctuating  samples  N 

3.  Chi  Square,  k=2;  Correlated  over  N  2 

4.  Chi  Square,  k*»2;  Fluctuating  over  N _ ZN. 


3.  CFA  -  2.36  t-log10  PFA]'<  -  1.02 
T  *  1.8  PD  -  0.9 
Gd  =■  1.23T/ (1  -  T2)*5 

x  “gfa+gd 

Y  -1/K 

LF  =  {-ln(PD}  G~  } 

FA 

x2l 

S  =  U  +  [1  +  16N/C5X2)]15} 

SDB  "  10  108lO(S) 


4.  No  analytic  inverse  for  Brooks'  formula  is  known. 


Table  A-2  Neuvy's  Algorithm  for  Minimum  Required 
Signal-to-Noise  Ratio 


Reference :  J.  Neuvy,  "An  Aspect  of  Determining  the  Range  of  Radar  Detection," 
IEEE  Trans.,  AES-9,  July  1970;  p  514-521. 


Definitions : 


*  Probability  detection  occurs,  given  that  target  signal 
is  present. 

■  Probability  of  false  alarm,  per  cell  per  N-sample  group. 
=  Effective  number  of  samples  integrated  at  video. 

■  Ratio  of  target  signal  energy. 


Algorithm: 


1.  E  -  exp{-N/3) 


2 .  Target  Distribution  Case _ a _ 

0.  Nonfluctuating  1  +  2E 

1.  Rayleigh;  Correlated  over  N  Samples  2/3(l+2/3E) 

2.  Rayleigh;  Uncorrelated  1 

3.  Chi-square,  k=2;  Correlated  over  N  Samples  3/4(l+2/3E) 

4.  Chi-square,  k32;  Uncorrelated _ 1 _ 


1/6+E 

2/3 

1/6+2/3E 


3-  S  -  a  log10(l/PFA)/{N2/3(logl0(l/PD))S} 


SDB  =  10  loSlO(S) 

,  _  rQ  loS10(1/PFA)  ,1/3 

4.  Inverse:  U  =  [ - jjz - ] 

S  N' ' 


PD  =  10 


APPENDIX  B 

Coefficients  of  Khoury-Bucknam 
Representations 


Definitions : 


i  Pn  ■  Probability  of  Detection 

I  D 

F.  SDB  ■  Minimum  Detectable  Signal- to-Noise  Ratio  in  dB 

L 

1  N  *  Number  of  signal  samples  integrated  at  video 

¥  P_.  =  Probability  of  False  Alarm 

i  rA 

Case™  Target  Distribution  Case 
'  0  ■  Nonfluctuating 

1  =  Rayleigh,  Correlated  across  N  samples 

2  =  Rayleigh,  Decorrelated  by  sample 

3  ®  Chi-square  with  k“2 ,  Correlated  across  N  samples 

4  »  Chi-square  with  k=2,  Decorrelated  by  sample 


Formulas : 


O.HS  +  AD) 

pd  -  i  -  (i  -  pfa)/[1  +  do  DB  )B]C 


„  10  ,  r,L"PFAa/C  ,  ! 

SDB  =  T  logio  i(T^}  -  1}  " 


i 

* 

I 

i 


Table  B-l  Coefficients  of  JChoury-Bucknam  Representations  Used 

in  Comparisons 


Distribution  Case 


10~6  10 


0.0433 

-13.64 

3.3177 

4.8476 

0.2938 

-5.32 

4.1934 

3.7524 


0.1265 

-8.98 

2.6832 

0.3205 

0.7345 

-1.34 

2.9222 

0.2962 


Case  1 


0.3936 


3.7306 

1.2269 


0.0912 

-10.40 

2.6786 

0.5805 

0.5309 

-2.75 

2.9538 

0.5250 


Case  3 


0.3475 

-4.59 

3.872 

1.8606 


10  100 


10~8  10 


1.3366 

+1.26 

5.0991 

3.3728 

0.0380 

-14.20 

4.5679 

2.7339 

0.263 

-5.80 

5.5040 

2.1039 


3.0761 

+4.88 

3.0927 

0.2815 

0.0991 

-10.04 

3.0893 

0.2742 

0.5861 

-2.32 

3.0900 

0.2778 


1.349 

+1.30 

4.9281 

3.0079 


Case  1 


0.3373 


4.2905 

1.0263 


2.4660 

+3.92 

3.4707 

0.4316 

0.0760 

-11.19 

3.2739 

0.4504 


0.4395 


3.2360 

0.4691 


1.3677 

+1.36 

5.0075 

3.0079 


Case  3 


0. 3090 
-5.10 
4.6284 
1.3621 


10  100 


10_L0  1 


1.0617 

+0.26 

6.3879 

3.3728 

0.0337 

-14.73 

4.6590 

3.0079 


2.917 
+4 , 64 
4.8264 
0.1733 

0.0796 

-10.99 

3.1868 

0.2670 


1.2359 

+0.92 

5.8621 

2.2816 


Case  1 


1.9953 

+3.00 

3.5950 

0.4221 

0.0590 

-12.29 

3.2712 

0.4598 


1.1749 

+0.70 

5.9854 

2.9219 


Case  3 


Table  B-l  (continued) 


Distribution  Case 


PFA 

N 

Coef . 

0 

— 

1 

2 

3 

4 

o 

t— 1 

1 

o 

H 

10 

A 

0.2317 

0.5164 

0.2897 

0.3724 

0.2786 

AD 

-6.35 

r— 

oo 

t 

-5.38 

-4.29 

-5.55 

B 

5.6485 

3.4370 

4.6576 

3.3817 

5.1942 

C 

2.6065 

0.2460 

0.9326 

0.4504 

1.1766 

IQ'10 

100 

A 

1.064 

2.2803 

1.1535 

1.7418 

1.1588 

AD 

+0.27 

+3.58 

+0.62 

+2.41 

+0.64 

B 

6.9134 

3.4561 

6.8119 

3.6463 

7.1206 

C 

2.8617 

0.2487 

1.7367 

0.4127 

1.8606 

APPENDIX  C 

Curve  Fitting  ?D  vs  SNR  Curves 

by  E.  N.  Khourv 
TSC  Memorandum  TSC-W21-36 


C-l 


TSC-W21-36/rad 
March  6,  1979 
A2095111 


Curve  Fitting  P^  va.  SNR  Curves 
by  E.  N.  Khoury 


The  probability  of  detection  (P  )  for  a  given  signal-to-noise 

d 

ratio  (S)  is  a  complex  function  of  the  target  mean  cross  section,  target 
fluctuation  statistics,  number  of  pulses  integrated  and  probability  of 
false  alarm  (detection  threshold  setting) .  Computing  this  function  re¬ 
quires  evaluating  some  rather  complicated  integrals  with  limits  which 
may  extend  to  infinity.  While  this  may  be  accomplished  by  using  a 
Simpson's  rule  integrator  or  other  such  techniques  in  an  iterative 
manner,  these  techniques  are  very  general  and  provide  accuracy  at  the 
expense  of  excessive  use  of  computer  time  and  core  storage  requirements. 

To  reduce  these  core  and  time  requirements  on  the  program,  an  approxi¬ 
mation  method  has  been  developed.  This  was  done  by  computing  a  fit  to  the  P^ 
vs.  SNR  case  desired  and  using  this  fit  in  the  Statistical  Report  Gener¬ 
ator.  The  curve  used  to  fit  the  data  points  is  a  three  parameter  low 
pass  filter  function  of  the  form 


vs) 


1  - 


_ 1 _ 

*  Skz  +  l]k3 


(1) 


where  P^(S)  is  the  estimate  of  ?d  for  any  given  value  of  SNR 

and  k^,  k^,  k^  are  the  parameters  which  are  varied  to  obtain  the 

best  fit  to  the  desired  model;  i.e.,  Swerling  Case  I,  M  pulses  integrated 

and  ?.  -  10-6. 

ta 

If  the  cost  function  is  selected  as  the  squared-error  criteria, 
the  function  to  be  minimized  is 


e2(Pd) 


N 

E 

i-1 


{pdi(V  - 


-  I  (P4<(S.)  -  [1  -  k  S  k*  +  l)~K5j}2 
i-1  1  11 


(2) 


1 


where 


fdicsi> 


is  che  number  of  data  points  selected  from  the  appro¬ 
priate  curve  of  vs.  SNR.  Not-'  that  these  points 
should  cover  the  entire  range  of  P^  values. 

is  the  true  value  of  for  selected  from  the 
appropriate  curve. 


P,.(S.)  is  che  estimate  of  P,.(S). 

Qi  1  Cl 

are  che  model  parameters. 


The  minimization  is  accomplished  in  a  standard  manner  by  taking 
the  partial  derivatives  of  the  squared-error  with  respect  to  the  three 
parameters  and  equating  to  zero.  These  partial  derivatives  are: 


3[e  (Pd). 
Sk, 


N 


*3  £  PdlSlk‘lkl  Slki  + 


+  2k,  E  S  k*  [k  S  kz  +  l)~k3-1 
J  i-1  1  1  1 


2  k,  Z  S  kj  (k  S  kz  +  l]"2  kj-1  -  C 
J  i-1  A  1 


(3) 


(*d> 

Sk„ 


-2  k,k,  Z  Pdi  lnCS^  S^2  [k^2  +  l]"’*3"1 


J  i-1 


+  2  kxk3  Z  In  (St)  S^2  (k^S^2  +  l]"k3~l 


-  2  kLk3  E  In  (St)  S^2  [k^2  +  l]"2  k3_1  - 


i-1 


(4) 


The  parameters  and  fits  obtained  for  five  cases  used  In  this  study 
are  presented  In  Tables  1  thru  5.  Note  Chat  the  column  labeled  "PDF" 

Is  the  fit. 
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APPENDIX  D 


Program  PDFI 


(0001 ) 
(0002) 
(0003) 
(0004) 
(OOOS) 
(0006) 
(0007) 
(0008) 
(0009) 
(0010) 
(0011  > 
(0012) 
(0013) 
(0014) 
(0015) 
(0016) 
(0017) 
(0018) 
(0019) 
(0020) 
(0021  ) 
(0022) 
(0023) 
(0024) 
(0025) 
(0026) 
10027) 
(0028) 
(00^9) 
(0030) 
(0031  ) 
(0032) 
(0033) 
(0034) 
(0035) 
(0036) 
(0037) 
(0038) 
(0039) 
(0040) 
(0041 ) 
(0042) 
(0043) 
(0044) 
(0045) 
(0046) 
(0047) 
(0048) 
(0049) 
(0050) 
(0051  ) 
(0052) 
(0053) 
(0054) 
(0055) 
(0056) 
(0057) 
(0058) 
(0059/ 

(006  j) 
(0061 ) 
(0062) 
(0063) 
(0064) 
(0065) 
(0066) 
(0067) 
(0068) 
(0069) 
(0070) 
(0071 > 
(0072) 


C  PROCRAM  PDFIT 

C 

c 

^insert  syscqm:>keys.  f 

c 

c 

IMPLICIT  DOUBLE  PRECISION  (A-H,  Q-Z) 

C 

INTEGER  FILE6(3),  CODE 

REAL  CASES, PFAS,  THS(3>,  SNRS(12),  PDS(12>,  PDFS(12>, 

1  ERRORS( 12) , ERRMS, SDERRS, SEC  1 , SEC2, SECS, SES 
INTEGER  I AR ( 15) 

C 

DIMENSION  PD( 12) ,  ED( 12), WK(3).  ISTOP( 12),  ST( 1 2) ,  ERROR ( 12 ) 
DIMENSION  SNR ( 1 2 ) , PDF ( 1 2 ) , TH < 3 ) ,  DEL (3),  DINC ( 3 ) , W ( 3, 12) 
DIMENSION  IPO),  STP(3),  NP(IOO) 

C 

COMMON /COM /SNR ,  PDF,  TH,  DEL,  DINC,  W,  M,  KS 
C 

DATA  WK/O.  5D0,  0.  5D0,  0.  5D0/ 

DATA  ST/O.  OOlDO,  0.  001D0,  0  001D0,  0.  OOIDO,  0.  OOIDO,  0.  001D0, 
10  OOIDO,  O.  OOIDO,  0.  OOIDO,  0.  OOIDO,  0.  OOIDO,  0.  OOIDO/ 

DATA  IPF/: 100600/ 

C 

C 

C  PRINT  INFORMATIVE  MESSAGES. 

C 

WRITE! 1,  3000) 

WRITE(1,  3100) 

C 

C  ENTER  INPUT  DATA. 

C 

WRITE! 1, 3200) 

READ ( 1 ,  3300 )FILE6 
WRITE! 1,  3400) 

READ( 1,  *) ICASE 

GO  TO  (1111,2222,3333),  ICA5E+1 
1111  WRITE! 1,  3500) 

READ ( 1 , * ) PFAS 
GO  TO  10 

2222  WRITEd,  3600) 

READd.  *>PFAS 
GO  TO  10 

3333  WRITEd,  3700) 

READd,  * ) PFAS 
10  WRITEd,  3800) 

READd  -  * )  NCASE 
WRITEd,  3900) 

READ ( 1, *) (NP< III >,  1 1 1  =  1, NCASE) 

WRITEd,  4000) 

C 

C  OPEN  OUTPUT  FILE. 

C 

CALL  SRCH**(K$WRIT, FILE6,  6,  2,  0,  CODE) 

IF (CODE.  NE.  0)  GO  TO  9999 
C 

C  INTEGER/SINGULAR  TH  REAL/DOUBLE  PRECISION  CONVERSION. 

C 

CASE=FLOAT( ICASE) 

PFA=DBLE(PFAS) 

C 

C  PROCESS  SWERLING  CASE  NCASE  TIMES. 

C 

DO  999  ICNT»1, NCASE 
C 

C  PRINT  CASE  NUMBER  BEING  PROCSSED. 

C 

WR I TE ( 1, 9000) ICNT 
C 

NNP=NP( ICNT) 

C 

C  INITIALIZATION. 


.ilBlLlLdixM^  -.111111  HI.. . wi  ' 


(0073) 
(0074) 
(0073) 
(0076) 
(0077) 
(0078) 
( 0079 ) 
(0080) 
(0081 ) 
(0082) 
(0083) 
(0084) 
(0083) 
(0086 ) 
(0087) 
(0088) 
(0089) 
(0090) 
(0091 ) 
(0092) 
(0093) 
(0094) 
(0093) 
(0096) 
(0097) 
( 0098 ) 
(0099) 
(0100) 
(0101 ) 
(0102) 
(0103) 
(0104) 
(0103) 
(0106) 
(010/) 
(0108) 
(0109) 
(0  10) 
(0111 ) 
(0112) 
(0113) 
(0114) 
(0113) 
(0116) 
(0117) 
(0118) 
(0119) 
(0120 ) 
(0121  ) 
(0122) 
(0123) 
(0124) 
(0123) 
(0126) 
(0127) 
(0129) 
(0129) 
(0130) 
(0131  ) 
(0132) 
(0133) 
(0134) 
(0133) 
(0136) 
(0137) 
<0138  ) 
(0139) 
(0140) 
(0141  > 
(0142) 
(0143) 
(0144  ) 


20 


C 

c 

c 

23 


L-3 

n-o 

IP  (  1 >-0 
KS-O 
IP (2>-0 
IP  (3>-0 
ISS-0 

DEL  (  1  >-0.  000 100 
DELC21-0.  000 IDO 
D£L(3>-0.  000 IDO 
STP  (  1  >-0.  OOOOOIDO 
STP (2>-0.  00000 1  DO 
STP ( 3) *0.  OOOOOIDO 
i sun-o 
on IN- 1 .  OD-2 
WEIGH- 1 .  ODO 
SK— 1 .  ODO 
DO  20  1-1.  12 
ISTOP ( I ) — 0 
CONTINUE 

SC ASE-DBUE ( FLOAT ( NVP ) ) 

IF  ( DABS  (CASE— l .  ODO).  LE.  0.  1  DO )  SCASE— 1 .  ODO 
IF  (CASE— 1 .  000.  LE.  -O.  IDO)  SCASE- 1 .  1D9 


COMPUTE  AOAPTIVE  INITIALIZATION  FOR  SWERLIN©  CASE  NUMBER  1. 
rn-dble ( float i nnp  >  > 

IP( CASE.  NE.  1  ODO)  00  TO  23 
TH  (  1  )  -0.  11 00*RN**0.  72DO 
Th  ( 2  )  -2.  12D0 
TH  (  3  )  — 0.  400 
CO  TO  27 

COMPUTE  ADAPTIVE  INITIALIZATION  FOR  SWERLIN©  CASE  NUMBER  2. 

CONTINUE 

TH(  1  >-0.  11D0*(DBL£<FL0AT(NNP)  )  )**0.  3 ADO 
DEVI  I—  DL0010(PFA)-6.  ODO 
SC— DSIQNd.  ODO.  DEVI  1) 

D2P-S©*(0.  0 1  DO ) *-* ( 1 .  0D0/RN)4(DAB8(DEV1 1 )  >•*<  1.  ODO-1.  0D0/RN1/2.  ODO 
D3P— — O. 8D0-D2F 

TH ( 2 ) - 1  979D0+3  ODO* ( ( RN- 1 .  ODO ) /RN ) **6.  ODO 
TH(2>-TH(2)»0.  8D0*D2P 

TH(  3  )  — O.  434D0+1.  36D0»(  (RN-1.  0D0>/RN>**6.  ODO 
IF(SO.  LT.  0.  0D0>  TH ( 3 ) -TH <  3  ) ♦D3P 


IDO) .  AND.  (NNP.  LE.  6 )  > TH ( 2 )  — TH ( 2 ) ♦  1 
100)  .  AND.  (NNP.  LE.  6)  )  TH<  3)  -  *W(  3  >  *1 


IDO)  .  AND.  (NNP.  LE.  6  >  )  TH  ( 1  )  -  i  M<  1  )  /  ..  3Dw 


IF(  (CASE.  LE. 

IF (  (CASE.  LE. 

IF(  (CASE.  LE. 

27  PD ( 1 ) —0  03D0 
PD ( 2  >  — 0.  IDO 
PD  ( 3  )  — O.  2D0 
PD  ( 4  )  — O.  3D0 
PD  (  3  >  — 0.  400 
PD  ( 6 )  — 0.  3D0 
PD ( 7 ) — 0  600 
PD ( 8 )  — 0.  7D0 
PD  ( 9 )  — 0.  8D0 
PD(  10>-0.  9D0 
FD<  1 1 )— 0.  93D0 
PD<  12 )  — 0.  99D0 
SECS-O, 

SES-O. 

c 

C  COMPUTE  SNR. 

C 

SNRMIN— 13.  ODO 
SNRMAX-33  000 
DO  100  1-1.  12 

T- ( 2.  O00*PD (  I )  —  1 .  ODO)  /  l.  1  IDO 
OD-1.  231D0*T /DSQRT ( 1  000-T*T> 

8-2.  36D0*DSQRT ( -DLOO 10( PPA ) >-l.  0200 


ODO 

8D0 


(  0 1 45  > 
(0146) 
(0147) 
(0148) 
(0149) 
(0150) 
(0151  ) 
(0152) 
(0153) 
(0154) 
(0155) 
(0156) 
(0157) 
(0158) 
(0159) 
(0160) 
(0161  ) 
(0162) 
(0163) 
(0164) 
(0165) 
(0166) 
(0167) 
(0168) 
(0169) 
(0170) 
(0171  ) 
(0172) 
(0173) 
(0174) 
(0175) 
(0176) 
(0177) 
(0178) 
(0179) 
(0180) 
(0181  ) 
(0182) 
(0183) 
(0194) 
(0185) 
(0186) 
(0187) 
(0188) 
(0189) 
(0190) 
(0191 ) 
(0192) 
(0193) 
(0194) 
(0195) 
(0196 ) 
(0197) 
(0198) 
(0199) 
(0200) 
(C201 ) 
(0202) 
(0203) 
(0204) 
(0205) 
( 0206 ) 
(0207) 
( 0208 ) 
( 0209 ) 
(0210) 
(0211) 
(0212) 
(0213) 
(0214) 
(0215) 
(0216) 


XQ- ( G+GD ) *+2 

SN=X0/4.  QD0/D3LE ( FLOAT ( NNP ) >*( 1.  ODO+DSQRT < 1  0D0+16. ODO* 
1  D3LE( FLOAT (NNP ) ) /X0> ) 

IF(CASE.  EQ.  0.  ODO)  SLF-l.ODO 
IF  (CASE.  EQ.  0.  ODO)  00  TO  40 
Fl-( 1. ODO+GD/G ) * ( -DLQG ( PD ( I ) ) > 

IF  (CASE.  EQ.  1.  ODO)  GO  TO  30 

SLF=1.  0D0/(F1 >**(1.  ODO/DBLE ( FLOAT (NNP ) ) ) 

GO  TO  40 

30  SLF= 1 . 0D0/F1 

40  SN=SN*SLF 

SN-10. ODO*DLOG10(SN) 

SNRMIN=SN-2.  ODO 
SNRMAX»3N+2. ODO 
DO  80  J»l,  10000 
SN*(SNRMIN+SNRMAX)/2.  ODO 
CALL  TIMDAT ( IAR,  15) 

SEC1=FL0AT( IAR (7) ) +FLOAT ( IAR<8) ) /FLOA T ( I AR ( 1 1 ) ) 

CALL  DETECT  <  NNP,  FFA,  SN,  SCASE,  1 .  0D-4,  PROB  > 

CALL  TIMDAT< IAR, 15) 

SEC2=FL0AT( IAR ( 7) ) +FLOAT ( I AR < 8 ) > /FLOAT ( IAR ( 1 1 ) ) 
SECS=SECS+  <  SEC2-SEC 1 > 

D=PROB-PD( I ) 

IF ( D ) 50, 60, 70 

50  IF(DABS(D>.  LE.  DMIN)  GO  TO  90 

SNRM IN=SN 
GO  TO  80 
60  GO  TO  90 

70  IF(DABS(D) . LE.  DMIN)  GO  TO  90 

SNRMAX=SN 
80  CONTINUE 

90  SNR( I >=10. 0DO**(O. 1D0*SN> 

PD( I ) =PRQB 

SES=SES+SECS /FLOAT ( J ) 

SNRMAX=35. ODO 
100  CONTINUE 

SLi*SES/FLOATC I ) 

TH( 1 )«1.  DO/TH( 1 ) 

IF ( (CASE.  GE  9DO  ) . AND  (CASE. LE.  1.  IDO) )  TH C 1 >=SNR(3> 
IF(CASE.  LE.  .  IDO)  TH(1)»SNR(9) 

C 

c 

C  COMPUTE  INITIAL  DECUPLING  MATRIX. 

call  vbar 

CALL  DCPLE 
DO  160  KS=1, 1000 
n=n+i 
call  vbar 
c 

C  COMPUTE  DECUPLING  MATRIX. 

C 

IF(M.  GT.  L)  CALL  DCPLE 
C 

C  FORM  ERRORS. 

C 

DO  110  1  =  1,  12 

£D( I >=(PDF( I ) -PD ( I ) ) *WE I GH 

IF ( DABS ( ED ( I) ).  LE.  ST ( I ) )IST0P(I)  =  1 

I SUM* I SUM+ 1  STOP ( I ) 

110  CONTINUE 
C 

C  COMPUTE  WEIGHTED  ERRORS, 

C 

DO  130  1  =  1,  3 
D I NC ( I >=0, CDO 
DO  120  K=l, 12 

DINC ( I >=DINC ( I )+W( I,  K)*ED(K> 

120  CONTINUE 
C 

C  COMPUTE  UPDATED  PARAMETER  ESTIMATE. 

C 

TH ( I ) =TH ( I ) -WK ( I ) *D INC ( I >/SK 


(0217) 
(0218  > 
(0219) 
( 0220  ) 
(0221  > 
(0222) 
(0223) 
(0224  ) 
(0225) 
( 0226 ) 
(0227  ) 
(0228) 
(0229) 
(0230) 
( 023 1  ) 
( 0232 ) 
(0233) 
(0234) 
(0235) 
(0236) 
(0237  ) 
( 0230  ) 
(0239  ) 
(0240) 
(0241  ) 
(0242) 
( 0243  ) 
(0244  ) 
(0245) 
(0246  ) 
(0247) 
(0243  ) 
(0249  ) 
(0250) 
(0251  ) 
(0252) 
(0253) 
(0254  ) 
(0255) 
(0256) 
(0257) 
(0258) 
(0259) 
(0260) 
(0261  ) 
(0262) 
(0263) 
(0264  ) 
(0265) 
(0266 ) 
(0267) 
(0268) 
(0269) 
(0270) 
(0271  ) 
(0272 ) 
(0273) 
(0274  ) 
(0275) 
(0276) 
(0277) 
(0278) 
(0279) 
( 0290 ) 
(0281  ) 
<  0282 ) 
(0283 ) 
(0284) 
(0235 ) 
(0286  ) 
<0287 ) 
(0288  ) 


S3*WK( I )*DINC  (  I  ) 

I F  (  DABS  ( SS  )  LE.  ST i*  ( I  )  )  IP  ( I  >  *1 
I3S=IS3+IP( I  ) 

C 

C  CHECK  IF  ESTIMATE  HAS  CONVERGED,  IF  YES.  SET  ISTOP(I) 

C  EQUAL  TO  1. 

C 

130  CONTINUE 
C 

C  HAVE  ALL  THREE  ESTIMATES  CONVERGED? 

IF (  I  SUM.  EO.  12)  GO  TO  170 
IF(  ISS.  EQ.  3)  GO  TO  170 
DO  140  1=1. .2 
I STOP ( I )=0 
140  CONTINUE 

DO  1 50  1-1.3 
IP ( I )-0 

150  CONTINUE 
I5S-0 
I 5UM-0 

160  CONTINUE 
170  CONTINUE 
S2—0.  ODO 
SS-O.  ODO 
DO  180  1  =  1.  12 
ERROR ( I ) =PD ( I ) -PDF  < I ) 

S2-S2+ERR0R  < I >**2 
SS-SS+ERROR  C I ) 

180  CONTINUE 

ERRM-SS/12.  ODO 

3DERR=D3QRT ( S2/ 12.  ODO-ERRM*ERRM ) 

THlD-1  DO/TH  (  1  ) 

C 

C  DOUBLE  PRECISION  TO  SINGLE  PRECISION  CONVERSION. 

C 

CASES—SNGL ( CASE ) 

PFAS-SNGL(PFA) 

DO  190  1=1.3 

THS  < I ) —SNGL ( TH  < I )  ) 

190  CONTINUE 

DO  200  1  =  1,  12 

SNRS ( I ) —SNGL <  SNR  < I ) ) 

PDS ( I ) —SNGL ( PD ( I ) ) 

PDFS ( I ) —SNGL  <  PDF ( I ) ) 

ERRORS ( I ) -SNGL (ERROR ( I ) ) 

200  CONTINUE 

ERRMS=SNGL ( ERRM ) 

SDERRS-SNGL ( SDERR ) 

C 

C  WRITE  TO  OUTPUT  FILE. 

C 

WRITEC6, 1400) ICNT 
WRITE(6.  1420) 

WR I TE ( 6, 1440) 

WRITE (6. 1460) 

WR ITE(6, 1500) 

WRITE (6, 1600) 

WRITEC6. 1700)CASES, PFAS. NNP 
WR  ITE(6,  1800) 

WR  I  TE ( 6,  1900) 

WR I TE ( 6, 2000) 

WR  ITECq, 2100) 

WR  I  TE  (  6,  2200  )  THS  (  1  ) ,  THS  ( 2  ) ,  THS  *.3  ) 

WR ITE<6. 2300 ) 

WR  ITE(6,  2400) 

WR  ITE(6» 2500 ) 

WR  ITE(6. 2600) 

WR ITE(6, 2700) (I,  SNRS< I ) , PDS ( I > ,  PDFS( I ) ,  ERRORS ( I > .  1  =  1,  12) 
WR  ITE(6, 2800 ) ERRMS. SDERRS 
C 

C  GO  TO  TOP  OF  PAGE  AND  OBTAIN  NEXT  CASE. 


WR  I TE< 6,  2900) IPF 


<  0239 ) 
<0290) 

<  0291 ) 

<  0292 > 
(0293) 
<0294) 

<  0298 ) 
<0296) 
<0297) 

<  0290 ) 
<0299) 
<0300) 
<0301 1 
<0302) 
<0303) 
<0304) 
<0303) 
<0306) 
<0307) 
<0308) 
<0309) 
<0310) 
<031 1 ) 
<0312) 
<0313) 
<0314) 
<0313) 
<0316) 
<0317) 
(0318) 
(0319) 
(0320) 
<0321 ) 
<0322) 
<0323) 
( 0324 ) 

<  0323 ) 
<0326) 
<0327) 
( 0328 ) 

<  0329 ) 

<  0330 ) 
<0331 ) 

<  0332  I 

<  0333 ) 
! 0334 ) 
(0333) 
'0336) 
<0337) 

<  0338 ) 
(0339) 
(0340) 
<0341 ) 
(0342) 
(0343) 
(0344) 
(0343) 
(0346) 
(0347) 
(0348) 
(0349) 
(0330) 
<0351 ) 
(0332) 
(0353) 
(0354) 
(0333) 
(0336) 
(0337) 
(0338) 
<0339) 
<0360) 


CONTINUE 


WIND  OUTPUT  FILE. 

REWIND  6 

DSE  OUTPUT  FILE. 

CALL  SRCH*«(K»CLOS,  0.  0.  2,  0.  CODE) 
IF < CODE.  NE.  0 )  00  TO  9999 

00  TO  99999 


FORMAT  <  //2X,  'CASE  NUMBER'.  14) 

FORMAT  - - - -  - - - '  ) 

FORMAT <////36X,  'INPUT  PARAMETERS') 

FORMAT (36X.  ' - - '  ) 

FORMAT (  / /  1  3X.  'SWERLINO  CASE  NO.  '.  10X,  'PFA',  12X, 

1  'PULSES  INTEORATED') 

FORMAT <  1  3X.  ' - 10X.  ' 12X.  ' - ■ 

FORMAT ( //22X,  F3.  1,  1 1  X,  Ell.  1,  17X,  13) 

FORMAT < / // /34X.  'PARAMETER  ESTIMATES') 

FORMAT  <  34X . - — — - ') 

FORMAT(//10X.  'PARAMETER  l',  10X,  'PARAMETER  2'.  10X,  'PARAMETER  3') 

FORMAT (  1 8X .  ' - 10X,  ' - 10X,  ' - ') 

FORMAT!  /  /14X,  FI  4  7,  7X,  FI  4.  7,  7X,  FI  4.  7) 

FORMAT  (  /  /  /  /  42X,  'FIT') 

FORMAT ( 42X .  ' - ' ) 

F0RMAT(//19X,  'N',  10X,  'SNR',  10X,  'PD'.  10X,  'PDF',  10X,  'ERROR') 

FORMAT!  18X.  10X,  ' - ',  10X,  ' —  ',  1CX  ' - 10X,  - - - 

FORMAT!  17X.  12.  6X.  Ell.  4.  3X,  F8.  3.  4X,  F8.  3,  6X.  FS.  3) 

FORMAT  I///19X  'MEAN  ERROR  -  F7.  3,  12X,  'RmS  ERROR  ■  ' .  F7.  3) 

FORMAT !A2) 

F0RMAT<///7X.  'PROORAM  POFIT  COMPUTES  AN  F.'.'IRICAL  FIT  TO  THE'/, 

1  2X ,  '  EXACT  PD  VS.  SNR  CURVE  DESIRED.  THE  EXACT  CURVE  IS'/, 

1  2X.  '  SPEC  IFIEO  BY  SWERLINO  CASE  NUMBER.  PRDBA0LI TY  OF'/, 

1  2X,  '  FALSE  ALARM,  AND  NUMBER  OF  PULSES  NON-COHERENTLY ' /, 

1  2X,  '  INTEORATED.  THIS  PROORAM  WILL  RETURN  THE  THREE'/. 

1  2X.  '  PARAMETERS  WHICH  OIVE  THE  BEST  FIT  IN  THE  MMSE ' / , 

1  2X.  '  SEN8E  <CF  REF  TSC-W21-40)  ' /, 

1  6X ,  '  POFIT  WA8  DEVELOPED  BY  E.  N.  KHOURY,  AND  R.  F.  '/, 

1  2X , '  PUSATERI  -  JUNE  79  '///) 

FORMATOX,  'INPUT  DATA  CONSISTS  OF  THE  FOLLOW  I  NO:  '//, 


3X,  '  A 

3X,  '  0 

3X»  '  C 

3X,  '  D 

3X,  '  E 

FORMAT! //3X 
FORMAT ( 3A2 ) 
FORMAT' //3X 


A)  OUTPUT  FILE  NAME'/. 

0)  SWERLINO  CASE  NUMBER'/, 

C)  PRO0AALITY  OF  FALSE  ALARM'/, 

D)  NUMBER  OF  CASES  TO  BE  RUN'/, 

E)  NUMBER  OF  PULSES  INTEORATED  FOR  EACH  CASE') 
X,  'ENTER  OUTPUT  FILE  NAME  -  6') 


FORMAT' //3X,  'ENTER  SWERLINO  CASE  NUMBER  -  0  1  2') 

FORMAT! //3X, 'ENTER  PRQBABLITY  OF  FALSE  ALARM'/, 

1  '  RANOE:  1 .  E-3  THROUGH  l.E-13') 

FORMAT < / /3X,  'ENTER  PRO0ABLITY  OF  FALSE  ALARM  -  PFA'/, 

1  '  RANOE:  l.E-2  THROUOH  l.E-38't 

FORMAT! //3X, 'ENTER  PROBABLITY  OF  FALSE  ALARM  -  PFA'/, 
l  '  RANGE:  1 .  E-3  THROUOH  l.E-B'J 
FORMAT! //3X.  'ENTER  NUMBER  OF  CASES  TO  BE  RUN'/, 

1  '  RANOE:  1  THROUOH  100') 

FORMAT ( //3X,  'ENTER  NUMBER  OF  INTEORATED  PULSES  FOR  EACH  CASE'/, 
l  '  RANOE:  1  THROUOH  100') 

FORMAT!// ) 

FORMAT! 'PROCES9INO  CASE  NUMBER:  '-I4) 

FORMAT <  /2X i  023.  9,  3X,  D23.  8,  3X,  023.  8) 

FORMAT!  /2X,  12,  3X,  D23.  8) 

FORMAT (  /2X ,  13,  3X,  12,  3X,  D23.  S) 


<0361  ) 
(0362) 
(0363) 
(0364) 
(0363) 
<  0366 ) 


C 

9999  CALL  ERRPR* < K*NRTN, CODE, 0. 0. 


C 

C 

99999  CALL  EXIT 
END 


'PDFIT 


3) 


<0001  > 

SUBROUTINE  DCPLE 

<0002) 

C 

<0003) 

C 

<0004) 

IMPLICIT  DOUBLE  PRECISION  (A-H,  0-Z  > 

<0005) 

C 

<0006) 

DIMENSION  PS (12),  TS(  12),  A <3.  3)-  BC3,  3),  DEO,  12) 

<0007) 

DIMENSION  SNR ( 1 2 ) , PDF ( 1 2 ) ,  TH< 3 ) ,  DEL < 3 ) • DINC<3>, W(3,  12) 

<0008 ) 

c 

<000? ) 

COMMON/COM/SNR,  PDF, TH, DEL/  DINC* W» M#  II 

<0010) 

c 

<0011) 

c 

<0012) 

WEIGH=1.  DO 

<0013) 

DO  4  1-1,3 

<0014) 

DO  2  U=i, 3 

<00155 

A  <  I,  J >  *0. 

<0016) 

2 

CONTINUE 

<0017) 

4 

CONTINUE 

<0018) 

DO  10  1-1,  12 

<0019) 

PS< I >»PDF< I ) 

(0020) 

10 

CONTINUE 

(0021  ) 

DO  30  1  =  1,  3 

(0022) 

TS< I >=TH< I ) 

<0023) 

TH ( I ) *TH< I )+DEL< I  ) 

(0024) 

CALL  VBAR 

(0025) 

DO  20  J*l, 12 

(0026) 

DE  < I ,  J )  =  <  PDF ( J  > -PS ( J ) ) /DEL ( I ) /WE I GH 

(0027 ) 

20 

CONTINUE 

(0028) 

TH ( I ) =TS( I  ) 

<0029 ) 

30 

CONTINUE 

(0030) 

DO  60  MK-1,  3 

<0031  ) 

DO  50  1*1,  3 

(0032) 

DO  40  J=l,  12 

(0033) 

A  ( MK,  I  )  =DE  ( MK,  J  )  *DE  <  I ,  J  )  +  A  <  MK,  I  ) 

(0034) 

WC  I,  J )  =0. 

(0035) 

40 

CONTINUE 

(0036) 

50 

CONTINUE 

(0037) 

60 

CONTINUE 

(0038) 

CALL  MX  INVR  ( A,  B,  IBAD) 

<0039) 

IF(  IBAD.  EQ.  1  )  WRITE(  1,  1000) 

(0040) 

DO  90  MK*1,  3 

(0041  ) 

DO  80  J*l,  12 

(0042) 

DO  70  1*1,  3 

(0043) 

W(MK,  J ) *W(MK,  U)+B<MK,  I ) *DE ( I , J) 

( 0044 ) 

70 

CONTINUE 

<0045) 

00 

CONTINUE 

( 0046 ) 

90 

CONTINUE 

(0047) 

DO  100  1*1,  12 

(0040) 

PDF  < I )=PS< I ) 

( 0049) 

100 

CONTINUE 

(0050) 

GO  TO  108 

(0051  ) 

105 

CONTINUE 

<  0052) 

DO  107  1*1, 3 

(0053) 

DO  106  U=1 ,  12 

( 0034 ) 

106 

W<  I,  J ) =DE ( I, J) 

(0055) 

107 

CONTINUE 

(0056) 

108 

CONTINUE 

(0057) 

M«0 

(0058) 

C 

(0059) 

1000 

FORMAT<  'DCPLE:  **  ERROR  **  IBAD  *  1') 

(0060) 

C 

(0061 ) 

RETURN 

(0062) 

END 

<  OQOl  ) 
<0002)  C 

<  0003  >  C 

<  0004  ) 
(0005)  C 

<  0006  )  C 

<  0007  ) 

<  oooa  > 

<  0009 )  C 

<  0010  ) 
<0011  ) 


SU 3 ROUTINE  DETECT < N, PFA, SN, SWERL, TQL,  PROB ) 


IMPLICIT  DOUBLE  PRECISION  <A-H, 0-Z) 


T=THRESH(PFA, N, 32) 

PROB=PROBDE ( T < SN, SWERL, 32,  N. TOL> 

RETURN 

END 


<0001  > 

DOUBLE  PRECISION  FUNCTION  GFUNCT < T, N,  M ) 

<0002  ) 

C 

<  0003  > 

C 

<  0C04  ) 

C 

THIS  FUNCTION  COMPUTES  THE  INCOMPLETE 

<0005) 

C 

GAMMA  FUNCTION. 

(0006) 

C 

<0007) 

IMPLICIT  DOUBLE  PRECISION  <A-H, 0-2) 

<oooa> 

C 

<0009  ) 

COMMON /DETPD/  TO,  TT,  <5,  C, R.  D 

<0010) 

C 

<0011  ) 

C 

<0012) 

TO=T 

<0013) 

IF  (M.  GT.  20)  GO  TO  30 

(.0014) 

R»f1*N 

<0015) 

TT*T /DBLE<  FLOAT (M) ) 

<0016) 

C-l.  /<TT+1.  )**R 

<0017) 

D-TT/(1.  +TT ) 

<0018) 

G-0.  DO 

<0019 ) 

DO  20  K-l.  N 

( 0020  > 

G=G+C 

<0021  ) 

C  =  C  *  <  R+DDLE  <  FLOAT (K  > )-l.  ) *D/DBLE ( FLOAT < K ) > 

<0022) 

20 

CONTINUE 

<0023) 

GO  TO  50 

<0024) 

30 

TT=«T*DBLE<  FLOAT  <N>  > 

<0025) 

C=DEXP (-TT) 

<0026 ) 

G*0.  DC 

(0027) 

DO  40  K-l,  N 

<  0028 ) 

G—G+C 

<0029) 

C=C*TT /DELE  <  FLOAT (K) > 

<0030 ) 

40 

CONTINUE 

<0031  ) 

50 

GFUNCT=G 

<0032) 

C 

(0033) 

RETURN 

(0034  ) 

END 

(0001  ) 

1  SEQ 

NEST 

SUBROUTINE  MX INVR ( A.  8  IBAO > 

<  0003  > 

(0003) 

1 

SUBROUTINE  HXINVRtA.  B,  I8AD> 

(0004) 

(0005) 

3 

C 

SUBROUTINE  MXINVR  FACTORS  MATRIX  A  INTO  LU  CROUT  DECOMPOSITITON 

(0006) 

3 

c 

AND  COMPUTES  ITS  INVERSE.  MATRIX  B. 

(0007) 

(0000) 

4 

IMPLICIT  DOUBLE  PRECISION  (A-M.  0-Z> 

(0009) 

(0010) 

3 

c 

(001 1  ) 

(0013) 

6 

DIMENSION  AO.  3),  BO,  3> 

(0013) 

(0014) 

7 

c 

(0013) 

(0016) 

3 

FOR  <  t«l  TO  3) 

(0017) 

9 

1 

FOR  <  J*  1  TO  3) 

(0018) 

10 

2 

B(I,J)«0. 

(0019) 

1 1 

2 

IF  (  I.  EQ.  J)  8<  I.  J>»1. 

(0030) 

12 

1 

END  FOR 

(0031  ) 

13 

END  FOR 

( 0023 ) 

( 0033 ) 

14 

c 

(0024) 

(0033) 

13 

IBAD-0 

( 0026 ) 

16 

FOR (K«l  TO  3) 

(0037 ) 

17 

1 

FOR  (  I«K  TO  3) 

(0020) 

18 

3 

DINN1-0  DO 

(0039) 

19 

2 

F0R<L»1  TO  (K-l)> 

( 0030 ) 

30 

3 

DINN1-DINN1+A( I,L)*A<L, K> 

(0031 ) 

21 

3 

END  FOR 

(0032) 

22 

2 

A<  I,  K)-A(  I,  K1-DINN1 

(0033) 

23 

1 

END  FOR 

( 0034 ) 

( 0033 ) 

24 

1 

c 

*  TEST  FOR  3 1 NQLLAR  I TY. 

(0036) 

( 0037 ) 

23 

1 

IF(A<K,  K).  NE.  0) 

<  0030 ) 

26 

2 

FOR  < (K+l )  TO  3) 

( 0039 ) 

27 

3 

DINN3-0.  DO 

(0040) 

30 

3 

FORCL-l  TO  <K-1>> 

(0041 ) 

39 

4 

DINN2»DINN3*A<K,  L>*ACL. J) 

(0043) 

30 

3 

END  FOR 

(0043) 

31 

3 

A  (  K,  J)  ■  (  A  <  K,  J  )  -DINN2  )  /  A  (  K,  K) 

(0044) 

33 

2 

END  FOR 

(0043) 

33 

t 

ELSE 

( 0046 ) 

34 

2 

IBAD-l 

(0047) 

33 

2 

K»4 

( 0040  > 

36 

1 

END  IF 

( 0049 ) 

37 

END  FOR 

(0C30) 

(0031  ) 

38 

c 

(0032) 

(0033) 

39 

1FUBAD.  NE.  1  ) 

(0034) 

40 

1 

F0R1M-1  TO  3) 

(0033) 

41 

2 

FOR  < I " 1  TO  3> 

0036) 

43 

3 

DINN1-0.  DO 

(0037) 

43 

3 

FOR  ( IX"  1  TO  (I-l>> 

(0038) 

44 

4 

DINN1-DINN1*A( I,  K)*B(K,M) 

(0039) 

1  SEQ 

NEST 

SUBROUTINE  MXINVRIAi  B>  ISAD) 

i 0060 ) 

(0061 ) 

43 

3 

END  FOR 

(0062) 

46 

3 

B(  I,  M)«(B(  I.  M>-DINN1  ) /Al  I,  I) 

( 0063 ) 

47 

2 

END  FOR 

( 0064 ) 

48 

2 

FOR (1-3  TO  1  STEP  -1) 

(0063) 

49 

3 

DINN1-0.  DO 

( 0066  > 

30 

3 

FDR< J-( Ifl )  TO  3) 

(0067) 

31 

4 

DINN1«DINN1*A(  I,  U)*B(U.  M> 

(0068) 

32 

3 

END  FOR 

(0069) 

S3 

3 

8(  I.  M>-B<  I,  M)-DINN1 

( 0070 ) 

34 

2 

END  FOR 

(0071 ) 

33 

1 

END  FOR 

(0073) 

36 

ELSE 

<0073) 

57 

1 

RETURN 

(0074) 

58 

END  IF 

<0075) 

(0076) 

59 

c 

<0077  ) 
<0079) 

60 

RETURN 

<0079 ) 

61 

END 

<0030) 
<0031 ) 

IFTRAN 

STATI 

STICS 

<0082) 

61 

CARDS 

READ 

<0083) 

0 

ERROR (S)  FOUND 

(0001 ) 
<00021 

c 

(0003) 

c 

(0004) 

c 

(0003) 

c 

(0006) 

C 

(0007) 

c 

(0008) 

(0009) 

c 

(0010) 

(0011) 

c 

(0012) 

(0013) 

c 

(0014) 

(0013) 

c 

(0016) 

c 

(0017) 

(0018) 

(0019) 

c 

(0020) 
(0021  ) 

c 

( 0022 ) 

40 

(0023) 

(0024) 

(0023) 

(0026) 

( 0027 ) 

c 

(0028) 

c 

(0029) 

c 

(0030) 
(0031  ) 

( 0032 ) 
(0033) 

41 

( 0034 ) 

42 

(0033) 

(0036) 

(0037) 

30 

(0030) 

<  0039 ) 
(0040) 
(0041  ) 
(0042) 

c 

(0043) 

C 

(0044) 

c 

(0043) 

(0046) 

(0047) 

(0048) 

31 

(0049) 

32 

(0030) 
(0031  ) 
(0032) 
(0033) 

C 

(0034) 

(0033) 

C 

(0036) 

33 

(0037) 

(0038) 

r 

(0039) 

(0060) 

C 

(0061  ) 

c 

(0062) 

c 

(0063) 

60 

(0064) 

C 

(0063) 

( 0066 ) 

C 

' 0067 ) 

c 

(0068) 

10 

(0069) 

c 

(0070) 

DOUBLE  PRECISION  FUNCTION  PROBDE (T,  X09.  AH,  H,  N.  ERR ) 

FUNCTION  CALCULATES  PROBABILITY  OF  DETECTION  FOR 
OIVEN  THRESHOLD  T,  CA8E  AK,  NUMBER  OF  PULSESN. 

NUMBER  OF  CFAR  CELLS  M.  AND  SIONAL-T0-NO18E  XDB. 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-i) 

INTE8ER-4  j 

COMMON  /DETPD/  TO,  TT,  00,  CO,  R,  D 
DATA  T0/-1.0D0/ 

IF  (  DABS  I  AK-DBLE(FLOAT  (N )  )  )  OT.  0.  IDO )  00  TO  40 
PR03DE-0FUNCT (  T /  ( 1.  DO-MO.  DO**  < .  1D0*XDB)  ) .  N,  M) 

RETURN 

IF(T.  NF.  TO)  O-OFUNC T < T,  N,  M > 

O-CO 

C-CO 

XN-DBLEt FLOAT (N) )*10  DO** (  1D0*XDB) 

IFlAH.  LE.  1000.  DO)  00  TO  41 

CHECK  FOR  POSSIBLE  UNDERFLOW. 

IF(XN.  OT  700.  DO)  00  TO  60 
A-DEXP(-XN) 

00  TO  42 

A«(l  DO*XN/AK)«*(-AK) 

E-l  DO-A 
J-0 

PD»A»0 

A«A«XN/ ( 1 .  OO-DBLE  <  FLOAT <  J  >  >  > 

IF  <  AK.  LE.  1000.  00)  A»A*<1  DO*  ( DBLE ( FLOAT ( J  )  )  )  / AK  )  /  ( 1 .  DO+XN/AK) 

1 

IF<FLOAT(  J).  OT.  I.  E6)  00  TO  33 
0»C*C 

IF  M.  OT.  20.  TREAT  AS  IDEAL  CFAR. 

IF(M.  OT.  20)  00  TO  31 

S=C5  ‘  R*2Bt-E  < FLOAT  <  N* ( FLOAT  ( J  )  )  -1)  )  )  *0/DBLE  ( FLOAT  <  N*  ( FLOAT  ( J )  )  )  ) 
00  TO  32 

C -C*TT/ DBLE  <FL0AT<N4< FLOAT  C  J ) ) ) ) 

PD-PD*A*0 

E-E-A 

IF  <E.  OT  ERR  )  80  TO  30 
PROBDE-PD 

RETURN 

continue 

PROBDE-PD 

RETURN 

PRINT  ERROR  MESBA8E  -  UNDERFLOW. 

WRITE! 1,  1000) 

STOP 

>0  FORMAT! /2X.  33HUN0ERFL0W  FOR  NONFLUCTUATIN©  CASE) 

END 


<0001  ) 
<  0002  > 
<0003) 
( 0004 ) 
<  0005 ) 
<0006) 
<0007) 
<0008 ) 
<0009) 
<0010  ) 
<0011  ) 
(0012) 
(0013) 
(0014) 
(0015) 
(0016) 
(0017) 
(0018) 
(0019) 
(0020) 
(0021  ) 
(0022) 
(0023) 
( 0024  > 
( 0025 ) 
(0026) 


DOUBLE  PRECISION  FUNCTION  THRESH < PFA,  N.  M) 

C 

c 

C  THIS  function  computes  the  detection  threshold 

C  FOR  A  GIVEN  PROBABILITY  OF  FALSE  ALARM  PFA, 

C  NUMBER  OF  PULSES  N,  AND  NUMBER  OF  CFAR  CELLS  M. 
C 

IMPLICIT  DOUBLE  PRECISION  (A-H,  O-Z) 

C 

C 

ERR=PFA/100.  DO 
T0*1 . 

T  1  *1 .  1 

F0=GFUNCT<  TO, N, M)-PFA 
DO  20  1=1, 1000 
F 1 =GFUNC  T ( T 1 , N, M)-PFA 
IF(DABS(F1 ).  LT.  ERR)  GO  TO  30 
T2“T 1 -F 1  * ( T 1-TO  > / ( FI -FO ) 

TO=T  1 
T 1  =T2 
FO=F  1 

20  CONTINUE 

30  THRESH=T 1 

C 

RETURN 

END 


<0001  ) 
<0002)  C 
(0003)  C 
<0004 ) 
<0005)  C 
<0006) 
<0007)  C 
<0008  ) 

<  0009  )  C 
<0010)  C 
<0011) 

<0012 ) 
<0013) 

<0014  ) 

<001 5) 
<0016)  10 
<0017)  C 
<0013) 

<0019  ) 


SUBROUTINE  VBAR 


IMPLICIT  DOUBLE  PRECISION  <A-H, 0-2) 

DIMENSION  SNR  <  1 2 ) ,  P DF ( 1 2 ) i TH <  3  > , DEL  < 3  > ,  DINC  <  3 ) ,  W < 3,  12) 
COMMON/COM/SNR,  PDF,  TH,  DEL,  DINC,  U.  M.  1 1 


DO  10  1=1, 12 

IFIPDfTu ‘  i :  i ! :?N& 1  -  ’ ~™< 2 > 

CONTINUE 


RETURN 

END 


!  0001  ) 

C 

(0002) 

(0003) 

(0004) 

c 

(OOOS) 

(0006) 

c 

(0007) 

c 

! 0008 ) 

c 

(0009) 

(0010) 

c 

(0011  > 

c 

(0012) 

c 

(0013) 

(0014) 

(0015) 

c 

(0016) 

30 

(0017) 

c 

(0013) 

c 

(0019) 
(0020) 
<0021  ) 
(0022) 
(0023) 
(0024) 

c 

(0025) 

( 0026 ) 
(0027) 

55 

(002S) 

(0029) 

60 

(0030) 

65 

(0031  ) 

c 

(0032) 

c 

(0033) 

(0034) 

(0033) 

C 

(0036) 

100 

(0037) 

C 

(0038) 

C 

(0039) 

(0040) 

c 

<0041 ) 

c 

(0042) 

c 

(0043) 

(0044) 

c 

(0045) 

c 

(0046) 

c 

(0047) 

(0043) 

(0049) 

(0050) 

c 

(0051  ) 

c 

(0052) 

c 

(0053) 

(0034) 

(0055) 

(0056) 

(0057) 

(0058) 

(0039) 

(0060) 

c 

(0061  ) 

c 

(0062) 

c 

(0063) 
(0064) 
(0065 ) 
(0066) 

c 

(0067 ) 

300 

(0068 ) 
(0069) 

c 

(0070) 

1001 

(0071 > 

1002 

(0072 ) 

1003 

PROGRAM  QFIT 

DIMENSION  S(50),  SDB (30),  POO).  PFI T ( 30 ) .  ERR  <  SO).  RELERR  (  SO ) 
DIMENSION  SDBW  (  IS),  WOO),  AK  ( 3  ) .  PDW  <  50  > ,  TH<3> 

PF<A1.  B 1 .  C 1  >  X  >  *1 .  —  ( 1 .  -M  A1*X )**B1  )**(-Cl  ) 

GET  INPUT  DATA. 

CALL  QFIN(N, SDB , P) 

CONVERT  FROM  DB  TO  POWER  RATIO. 

DO  SO  1=1. N 

S  (  I  )  -10.  **  (  SDB  (  I  >  *  1  > 

CONTINUE 

ESTIMATE  COEFFICIENTS. 

WRITE!  1,  1002) 

READ ( 1. *> I USE 
IFUUSE  .  NE.  1  )G0  TO  55 
CALL  ESTCOF < N, P , SDB , AK  > 

GO  TO  63 
DO  60  1*1,  N 

PDW(  I  )  *=P  <  I  >*100. 

SDBW ( I )=SDB < I ) 

CONTINUE 

CALL  QPDFIT<PDW,  SDBW,  N,  AK > 

CONTINUE 

COMPUTE  FITED  VALUES  OF  PD. 

DO  100  1*1,  N 

PFIT(  I  )=»PF  <  AK(  1  )  ,  AK  (  2  )  .  AK  (  3  >  ■  S  <  I  >  ) 

CONTINUE 

COMPUTE  ERRORS. 

CALL  ERRORS (P, PFIT,  N, ERR,  RELERR, AVERR, AVREL ) 

OUTPUT  PRINTED  RESULTS. 

CALL  QFQUT ( SDB,  S,  P.  PFIT,  ERR,  N,  AK) 

PLOT  RESULTS  IF  DESIRED. 

WRITE! 1, 1001 ) 

READ! 1, *) IPLOT 

IF!  I  PLOT  .  EQ.  1  )  CALL  FTPLOT  ( N,  S,  P,  AK ) 

DO  COMPLETE  PDFIT  IF  DESIRED, 

WRITE! 1, 1003) 

READ! 1, *> IPDF 
IF! IPDF  .  NE.  1 )G0  TO  300 
TH! 1 )=AK( 1 ) 

TH ( 2  >  =AK ( 2 ) 

TH ( 3  >  =AK ( 3 ) 

CALL  PDFIT!N, SDB,  P,  TH) 

PLOT  RESULTS  IF  DESIRED. 

WRITE! 1.  1001 ) 

READ!  1,  *>  IPLOT 

IF  !  IPLOT  .  EQ.  1  ) CALL  FTPLOT !  N,  S,  P,  TH  ) 

CONTINUE 

CALL  EXIT 

FORMAT!  'IF  A  PLOT  OF  RESULTS  IS  DESIRED  ENTER  1,  ELSE  O') 
FORMAT!  -TO  USE  BOBS  VERSION  ENTER  1 ,  ELSE  O') 

FORMAT!  'TO  DO  COMPLETE  PDFIT  ENTER  1,  ELSE  O') 


(0073) 


END 
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000000 
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000000 
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000000 

0006 
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0034 
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"1002 

002103 

0020 

007  ID 

1003 

002127 

0054 
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“300 

002047 

0056 

0067D 

“50 

001577 
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55 
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“60 

001651 
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(0001  ) 

SUBROUTINE  DCPLE 

(0002) 

c 

(0003) 

c 

(0004) 

IMPLICIT  DOUBLE  PRECISION  (A-H, 0-2) 

(000S) 

c 

(0006) 

DIMENSION  PS<  12)  ,  TS  (  12),  A  (3,  3),  B(3,  3),  DEO,  12) 

(0007) 

DIMENSION  SNR ( 1 2 ) , PDF ( 1 2 ) , TH(3> . DEL ( 3 ) , DINC  (3) ,  W<3. 12) 

(0008) 

c 

( 0009 ) 

COMMON/COM/SNR,  PDF,  TH,  DEL, DINC,  W,  M.  II-  JV 

(0010) 

c 

(0011  ) 

c 

(0012) 

WE  I  OH*  1.  DO 

(0013) 

DO  4  1*1,  3 

(0014) 

DO  2  J*l,  3 

(0015) 

A (  I.  J )  =0. 

(0016) 

2 

CONTINUE 

(0017) 

4 

CONTINUE 

(0018) 

DO  10  1=1, JV 

>0019) 

PS( I ) =PDF ( I ) 

( 0020 ) 

10 

CONTINUE 

(0021  ) 

DO  30  1  =  1,3 

( 0022 ) 

TS ( I )=TH( I  ) 

( 0023 ) 

TH ( I ) =TH  < I ) +DEL (  I  ) 

( 0024 ) 

CALL  V3AR 

(0025) 

DO  20  J=l, JV 

( 0026 ) 

DE( I, J)=(PDF( J)-PS( J) ) /DEL ( I > /WEIGH 

( 0027 ) 

20 

CONTINUE 

( 0028 ) 

TH ( I )=TS< I ) 

(0029) 

30 

CONTINUE 

( 0030 ) 

DO  60  MK=1, 3 

(0031  ) 

DO  50  1  =  1,3 

( 0032 ) 

DO  40  J=l, JV 

(0033) 

A ( MK« I )=DE(MK, J ) *DE ( I , J)+A(MK, I ) 

(0034) 

W(  I,  J  )  =0. 

(0035) 

40 

CONTINUE 

(0036) 

50 

CONTINUE 

(0037) 

60 

CONTINUE 

( 0038 ) 

CALL  MX INVR (A, B, IBAD) 

( 0039 ) 

IF(  IBAD.  EQ.  1  )  WRITE(  1,  1000) 

(0040) 

DO  90  MK=1, 3 

(0041  ) 

DO  80  J=l,  JV 

(0042) 

DO  70  1  =  1,  3 

( 0043 ) 

W(MK, J)=W(MK, J)+B(MK, I )*DE< I, J> 

(0044) 

70 

CONTINUE 

(0045) 

80 

CONTINUE 

( 0046 ) 

90 

CONTINUE 

(0047) 

DO  100  1=1, JV 

(0048) 

PDF ( I >=PS( I > 

(0049) 

100 

CONTINUE 

(0050) 

GO  TO  108 

(0051  ) 

105 

CONTINUE 

(0052) 

DO  107  1  =  1,3 

(0053) 

DO  106  J=l, JV 

(0054) 

106 

W( I, J ) =DE ( I, J) 

(0055) 

107 

CONTINUE 

(0056) 

108 

CONTINUE 

(0057) 

M=0 

(0058) 

C 

(0059) 

1000 

FORMAT (  'DCPLE:  **  ERROR  **  IBAD  =  1  •  ) 

( 0060 ) 

C 

(0061  ) 

RETURN 

(0062) 

END 

**■  "ffftlff . 
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<0001 ) 
(0002) 

<  0003 ) 
<0004) 
<000S) 
<0006) 
<0007) 
<0003) 
<0009) 
(0010) 
<001 1 ) 
(0012) 
<0013) 
<0014) 
<0019) 
<0016) 
<0017) 
(0019) 
<0019) 
(0020) 
(0021) 
(0022) 

<  0023 ) 
<0024) 
(0029) 
<0026) 
<0027) 
( 0028 ) 

<  0029 ) 

<  0030 ) 
<0031 ) 
<0032) 

<  0033 ) 

<  0034 ) 

<  0039 ) 
(0036) 

<  0037 ) 

<  0030 ) 

<  0039 ) 
(0040) 
<0041 ) 

<  0042 ) 

<  UW43 ) 

<  0044 ) 
(0049) 
<0046) 

<  0047 ) 

<  0040 ) 
<0049) 
<0090) 
<0091 ) 

<  0092 ) 

<  0093 ) 
<0094) 
<0099) 
(0096) 
(0097) 
(0098) 
(0099) 

<  0060 ) 
<  006 1  ) 
(0062) 
(0063) 
(0064) 
(0069) 
( 0066 ) 
(0067) 
(0068) 
(0069) 
(0070) 
<0071 > 
<0072) 

(0073) 

(0074) 

(0079) 


SUBROUTINE  ERRORSIXACT,  XHES.  N,  ERR.  RELERR.  AVERR.  AVREL.  BDVERR) 
DESCRIPTION 

COMPUTES  THE  ABSOLUTE  AND  RELATIVE  ERRORS  FOR  EACH  DATA  POINT. 
ALSO  COMPUTES  THE  AVER ACE  A 8 801 UTE  ERROR.  AVERAGE  RELATIVE  ERROR. 
Si  THE  STANDARD  DEVIATION  OF  THE  DATA. 

METHOD 

USES  THE  STANDARD  DEFINITIONS  FOR  ABSOLUTE.  RELATIVE  ERR0R8. 

It  STANDARD  DEVIATION. 

VARIABLE  DESCRIPTIONS 


NAME  TYPE  DIM 


AVERR  R 
AVREL  R 
ERR  R 
RELERR  P 
SDVERR  R 


DESCRIPTION 


NUMBER  OF  DATA  POINTS 
ACTUAL  VALUES  OF  DATA 
MEASURED  VALUES  OF  DATA 


AVERAOE  ABSOLUTE  ERROR 
AVERAGE  RELATIVE  ERROR 
ABSOLUTE  ERROR  OF  EACH  POINT 
RELATIVE  ERROR  OF  EACH  POINT 
STANDARD  DEVIATION  OF  DATA 


UNITS  VIA 


CALL  SEC 
CALL  SEO 
-  CALL  SCO 


CALL  SEQ 
7.  CALL  SEQ 
CALL  SEQ 
7.  CALL  9EQ 
-  CALL  SEQ 


IMPORTANT  LOCAL  VARIABLES:  NONE 
EXTERNAL  REFERENCES 


REFERENCE  BOOHS.  REPORTS/MEMOB 


R.  E.  BLASE 

DIMENSION  XACT<  1  ),  XMES<  l  > .  ERR  <  1  ) ,  RELERR  <  1  ) 

COMPUTE  ERRORS  FOR  EACH  COMPONENT. 

51- O. 

52- 0 
SSQ-O. 

IBAD*0 

DO  90  1-1,  N 

ERR  < I > "(ACT  < I ) -XME9  < I ) 

IF  <  XACT  (  I  )  .  EQ.  0.  )  IBAD-IBAD+I 

I F <  XACT (  I  1  NE.  0  ) RELERR!  I  ! •< ERR < I ) /XACT i I  M«100. 
81-81+ERR< I ) 

S2“S2*RELERR  < I ) 

SSQ"S8Q*ERR< I ) *ERR  < I ) 

CONTINUE 

COMPUTE  AVERAOE  ERRORS-. 

AVERR-S1/N 

NUSE-N-IBAD 

AVREL-S2/NU3E 

COMPUTE  STANDARD  DEVIATION 

SDVSRR*SQRT  <  380/N-AVCRR»AVERR ) 

RETURN 

END 


AVERR 

R 

ARGUMENT 

000010 

000  IS 

0066M 

0072 

AVREL 

R 

ARGUMENT 

000011 

000  IS 

006SM 

ERR 

R 

ARGUMENT 

000006 

0001S 

0047S 

0056M 

0058 

0059 

0061 

ERRORS 

R 

000000 

000  IS 

I 

I 

000205 

0055M 

0056 

0057 

0058 

0059 

0060 

IQAD 

I 

000212 

0054M 

0057M 

0067 

N 

I 

ARGUMENT 

000005 

000  IS 

0055 

0066 

0067 

0072 

NUSE 

I 

000213 

0067M 

0068 

RELERR 

R 

ARGUMENT 

000007 

000  IS 

0047S 

0058M 

0060 

SI 

R 

000216 

0051M 

0059M 

0066 

S2 

R 

000220 

0092M 

0060M 

0068 

SDVERR 

R 

ARGUMENT 

000012 

0001S 

0072M 

SORT 

R 

EXTERNAL 

000000 

0072 

SSQ 

R 

ARGUMENT 

000222 

0053M 

006 1M 

0072 

XACT 

R 

000003 

000  IS 

0047S 

0056 

0057 

0058 

XMES 

R 

ARGUMENT 

000004 

000  IS 

0047S 

0056 

_50 

000132 

0055 

0062D 

0000  ERRORS 
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<0001)  SUBROUTINE.  ESTCOFIN,  P,  SDB,  AK) 

<0002)  C 

<0003)  C  DESCRIPTION 

<0004)  C  - 

<0005)  C  COMPUTES  THE  1ST  ESTIMATE  OF  THE  PARAMETERS  FOR  THE 

<0006)  C  PDFIT  FUNCTION. 

<0007)  C 

<0008)  C  METHOD 

<0009)  C  - 

<0010)  C  USES  THE  ALOGRITHM  DEVELOPED  BY  J.  BUCKNAM  DESCRIBED  IN 

(0011)  C  A  MEMO  ON  'QUICK  PDFIT'. 

(0012)  C 

(0013)  C  VARIABLE  DESCRIPTIONS 

(0014)  C  - - 

(0015)  C 

<0016)  C  NAME  TYPE  DIM  DESCRIPTION  UNITS  VIA 

(0017)  C 

(0018)  C  INPUTS: 

<0019)  C  N  I  1  NUMBER  GF  <S,P)  DATA  PAIRS  -  CALL  SEQ 

(0020)  C  P  R  N  ACTUAL  PD  VALUES  -  CALL  SEQ 

(0021)  C  SDB  R  N  ACTUAL  SIGNAL-TO-NOISE  RATIOS  DB  CALL  SEQ 

<  0022 )  C 

(0023)  C  OUTPUTS: 

(0024)  C  AK  R  3  INITIAL  ESTIMATES  FOR  K1.K2,  K3  -  CALL  SEQ 

(0025)  C 

(0026)  C  LOCAL  VARIABLES: 

(0027)  C  A  R  1  SLOPE  OF  U  VERS  SDB  FIT 

(0028)  C  B  R  1  Y  INTERCEPT  OF  W  VESR  SDB  FIT 

<0029)  C  AK3UN  R  1  UNCORRECTED  ESTIMATE  FOR  K3 

<0030)  C  SRL  R  1  RATIO  OF  SLOPE  OF  TRUE  FUNCTION  TO  SLOPE 

(0031)  C  OF  ASYMPTOTIC  FUNCTION  AT  LARGE  PD 

(0032)  C  SRS  R  1  RATIO  OF  SLOPE  OF  TRUE  FUNCTION  TO  SLOPE 

<0033)  C  OF  ASYMPTOTIC  FUNCTION  AT  SMALL  PD 

(0034)  C  Y 1  R  1  1ST  VALUE  USED  FOR  SMALL  PD 

( 0035 )  C  Y2  R  1  2ND  VALUE  USED  FOR  SMALL  PD 

(0036)  C  YM  R  1  SLOPE  OF  Y  VERS.  SDB  CURVE 

<0037)  C  Wl  R  1  VALUES  OF  BUCKNAM  FUNCTION 

(0038)  C  Z3  R  1  1ST  VALUE  USED  FOR  LAROE  PD 

(0039)  C  Z4  R  1  2ND  VALUE  USED  FOR  LARGE  PD 

<0040)  C  ZM  R  1  SLOPE  OF  Z  VERS.  SDB  CURVE 

(0041)  C 

<0042)  C  EXTERNAL  REFS: 

(0043)  C  - 

<0044)  C  INSORT  -INSERTION  SORTING  ROUTINE 

<0043)  C 

<0046)  C  REFERENCE  BOOKS,  REPORTS/MEMOS 

<  0047 )  C  - 

<0048)  C  TSC-W7-73/RAD,  U.  BUCKNAM, 1979 

(0049)  C 

(0030)  C  AUTHOR 

(0031)  C  - 

<0032)  C  R.  E.  BLASE 

(0033)  C 

(0034)  DIMENSION  P ( 1 ) , SDB ( 1 ) , AK ( 3 ) 

(0033)  C 

<0036)  C  SORT  P  fy  S  ARRAYS  ACCORDING  TO  INCREASING  PD. 

<0037)  C 

(0038)  CALL  INSORT ( SDB, P,  1 ,  N ) 

(0039)  C 

(0060)  C  COMPUTE  SLOPE  FOR  SMALL  PD. 

<0061)  C 

<  0062 )  Y 1 -ALOG 1 0 ( -ALOG 10(1,  -P (  1  )  >  ) 

(0063)  Y2«AL0G 10 ( -ALOG 10(1.  -P<2>  )  ) 

(0064)  YM-(Y2-Y1 ) / ( SDB < 2 ) -SDB ( 1 ) > 

(0063)  C 

(0066)  C  COMPUTE  SLOPE  FOR  LARGE  PD. 

(0067)  C 

(0068)  NM1»N-1 

(0069)  Z 3= -ALOG 10 ( 1 .  -P  < NM1 ) ) 

(0070)  Z4*»-AL0G10(  1 .  -P  <  N)  ) 

(0071 >  ZM-( Z4-Z3)/(SDB(N)-SDB(NM1 ) ) 

(0072)  C 


<0073)  C 
<007'’,)  C 
<0075) 

<0076) 

<0077)  C 
<0078)  C 

<  007V )  C 
<0000) 

<0001  ) 

<  0002 ) 
<0083) 
<0084) 
<0085) 

<  0086 ) 

<  0087 ) 

<  0088 ) 

<  ^089 ) 

<  0090 ) 

<0091  ) 

<  0092 ) 
<0093) 
<0094) 

<0095)  C 
<0096)  C 
<0097)  C 
<0098) 

<  0099 ) 
<0100)  C 
<0101)  C 
<0102)  C 
<0103)  C 
<0104) 
<0105) 
<0106) 
<0107) 
<0108) 
<0109) 
<0110) 
<0111) 
<0112) 
<0113) 
<0114) 
<0115) 
<0116)  200 
<0117) 
<0118) 
<0119) 
<0120)  C 
<0121)  C 
<0122)  C 
<0123) 
<0124) 
<0125) 
<0126) 
<0127)  C 
<0128) 
<0129) 


COMPUTE  UNCORRECTED  VALUE  OF  K3. 

AK3UN=ZM/YM 
RAK3UN=1.  /AK3UN 

COMPUTE  CORRECTION  FACTOR. 

PTIL1-1  -< 1. -P< 1 ) )**<RAK3UN> 

PTIL2=1.  -  <  1.  -PC2)  ) ■** <  RAK3UN  ) 

PTIL3-1 .  -  < 1.  -P  <NM1 ) ) #*  < RAK3UN ) 

PTIL4*1.  -<  1.  -P<N>  ) **  ( R AK3UN  ) 

SRLT2-AL0G10<  <  <  1.  -PTIL4 )  *■*(  -RAK3UN )  >-l.  ) 
SRLT1»AL0G10<  <  <1.  -PTI L3 ) ** ( -RAK3UN ) )-l.  > 

SRLB2* < -RAK3UN ) *AL0G1 0 < 1 , -PTIL4) 
SRLB1*<-RAK3UN)+AL0G10< 1.  -PTIL3) 

SRL*  <  SRLT2-SRLT 1 > / ( SRLB2-SRLB 1 ) 

SRST2=AL0G 10  < <  ( 1,  -PT IL2 ) ** ( -RAK3UN ) )-l  ) 
SRST1*AL0G10C  <  < 1.  -PTIL1 ) ** ( -RAK3UN ) >-l.  ) 
CONST-ALOG( 10.  O) 

SRSB2-AL0G10C  < -CONST*R AK3UN > * < AL0G1 0 < 1 .  -PTIL2) ) ) 
SRSB1-AL0G10C  < -CONST*R AK3UN > * < ALOG 1 0 < 1 .  -PTIL1 ) ) ) 
SRS* <  SRST2-SRST 1 > / <SRSB2-SRSB1 ) 

CORRECT  K3. 

AK3*AK3UN*SRL/SRS 
RAK3=1.  /AK3 

COMPUTE  VALUES  OF  BUCKNAM  FUNCTION  FOR  GIVEN  PD'S. 
FIT  W  TO  W=A*S  +B 

SUM=FLOAT  <  N ) 

SUMX»0. 

SUMY*0. 

SUMX2-0. 

SUMXY=0 
DO  200  1*1, N 

SDB I *SDB  < I ) 

WI=AL0G10<  <  < 1.  -P< I > ) < -RAK3 ) )-l.  ) 
SUMX*SUMX+SDBI 
SUMY=SUMY+WI 
SUMX2*SUMX2+SDBI*SDBI 
SUMXY»SUMXY+SDBI*WI 
CONTINUE 

D-SUM*SUMX2-SUMX*SUMX 
B*<SUMX2*SUMY-3UMX*SUMXY) /D 
A=<SUMXY*SUM-SUMX*SUMY> /D 

COMPUTE  Kl,  K2 

AK10-10  *A 

AK  (  1  )  =>10.  **<B/AK10) 

AK  <  2  )  *»AK  1 0 
AK  <  3  )  =>AK3 


RETURN 

END 


A 

R 

000747 

0119N 

0123 

AK 

R 

ARGUMENT 

000006 

000  IS 

00943 

0124H 

0129M 

0126M 

AK10 

R 

000791 

0123M 

0124 

0129 

AK3 

R 

000793 

009SM 

0099 

0126 

AK3UN 

R 

000799 

0079M 

0076 

0098 

ALOO 

R 

EXTERNAL 

000000 

0091 

AL0010 

R 

EXTERNAL 

000000 

0062 

0063 

0069 

0070 

0084 

0089 

0006 

0087 

0089 

0090 

0092 

0093 

0111 

B 

R 

000797 

0118M 

0124 

CONST 

R 

000761 

009 1M 

0092 

0093 

D 

R 

000763 

0117H 

01 18 

01  19 

ESTCOF 

R 

000000 

000  IS 

FLOAT 

R 

EXTERNAL 

000000 

0104 

I 

I 

000769 

0109M 

01 10 

0111 

INSORT 

I 

EXTERNAL 

000000 

0093 

N 

I 

ARGUMENT 

000003 

OOOIS 

0109 

00 98A 

0068 

0070 

0071 

0003 

0104 

NM1 

I 

000770 

006SM 

0069 

0071 

0082 

P 

R 

arqument 

000004 

OOOIS 

00948 

00  98 A 

0062 

0063 

0069 

0070 

OOSO 

0081 

0082 

0083 

01  11 

ptili 

R 

000771 

0080M 

0090 

0093 

PTIL2 

R 

000773 

0001N 

0089 

0092 

PTIL3 

R 

000779 

0082M 

0089 

0087 

PT1L4 

R 

000777 

0083M 

0084 

0086 

RAK3 

R 

001011 

0099M 

0111 

RAK3UN 

R 

001013 

0076M 

0080 

0081 

0082 

0083 

0094 

0099 

0086 

0087 

0089 

0090 

0092 

0093 

30B 

R 

ARGUMENT 

000009 

OOOIS 

00948 

0098A 

0064 

0071 

0110 

SD0I 

R 

001019 

01 10M 

0112 

0114 

0119 

SRL 

R 

001017 

0088M 

0098 

SRLB1 

R 

001021 

0087M 

0088 

SRLB2 

R 

001023 

0086M 

0088 

SRLT1 

R 

001029 

0089M 

0088 

SRLT2 

R 

001027 

0084M 

0088 

SRS 

R 

001031 

0094M 

0090 

SRSBI 

R 

001033 

0093M 

0094 

6RS32 

R 

001039 

0092M 

0094 

SRST1 

R 

001037 

0090M 

0094 

SRSTS 

R 

001041 

0089H 

0094 

SUM 

R 

001043 

0 1  QAM 

0117 

0119 

SUMX 

R 

001049 

0109M 

01 12M 

0117 

0118 

0119 

SUMX2 

R 

001047 

0107M 

Oil  4M 

01  17 

0118 

SUMXY 

R 

001091 

0108M 

0119M 

0118 

0119 

SUMY 

R 

001093 

0106M 

0113M 

01  18 

0119 

UI 

R 

001099 

0111M 

0113 

0119 

VI 

R 

001097 

0062M 

0064 

Y2 

R 

001061 

0063M 

0064 

YM 

R 

001063 

0064H 

0079 

Z3 

R 

001069 

0069M 

0071 

Z4 

R 

001067 

0070M 

0071 

ZM 

R 

001071 

007 1M 

0079 

_200 

000627 

0109 

01160 
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(0001 ) 
(0002)  C 
(0003) 

( 0004 )  C 
(0005)  C 
<  0006 )  C 
(0007) 

( 0008 )  C 
(0009)  C 
(0010)  C 
(OOll ) 
(0012)  C 
(0013)  C 
(0014)  C 
(0015) 
(0016)  C 
(0017)  C 
(0018)  C 
(0019) 

( 0020 )  C 
(0021 ) 
(0022) 


SUBROUTINE  FTPLOT(N,  S,  P,  AK> 

DIMENSION  S(1 ) »  P  <  1  ) ,  AK  ( 3  > 

OPEN  PLOTER. 

CALL  PLOTS (4,  0, 0) 

PLOT  DATA. 

CALL  PTDATA(N<  S<  P ) 

PLOT  FIT 

CALL  PTQFIT(AK) 

CLOSE  PLOTTER. 

CALL  PLOT ( 0.  ,  0.  ,  999) 

RETURN 

END 


AK 

R 

ARGUMENT 

000006 

000  IS 

0003S 

0015A 

FTPLOT 

R 

000000 

000  IS 

N 

I 

ARGUMENT 

000003 

000  IS 

001 1A 

P 

R 

ARGUMENT 

OOOOOS 

000  IS 

0003S 

001 1A 

PLOT 

R 

EXTERNAL 

000000 

0019 

PLOTS 

R 

EXTERNAL 

000000 

0007 

PTDATA 

R 

EXTERNAL 

oooooo 

001 1 

PTQFIT 

R 

EXTERNAL 

000000 

0015 

S 

R 

ARGUMENT 

000004 

000  IS 

0003S 

001 1A 
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SEQ  NEST 
1 

2 

3 

4 

5 

6 

7 

8 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


C 

C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


50 


SUBROUTINE  INSORT(Z,W,  IF,  ID 
SUBROUTINE  INSORT ( Z,  W,  IF,  IL> 


<  Z, W>  ARE  DATA  PAIRS 

SORTING  W  TO  BE  IN  INCREASING  ORDER 

INSERTION  SORTING  USED 

IF  IS  THE  FIRST  ELEMENT  OF  ARRAY  TO  BE  SORTED 
IL  IS  THE  LAST  ELEMENT  OF  ARRAY  TO  BE  SORTED 


DIMENSION  Z< 1 >. W< 1 > 


USING  COMPILER  TIME  DOMAIN 
COMPUTATION  DONE  IN  PLACE 


ISTAR=IF+1 
DO  50  U-ISTAR, IL 


SUBROUTINE  TO  INSERT  W(IF+1>  TO  W(IL) 


CALL  INSERTC  Z,  W,  U,  IF,  IL) 
CONTINUE 
RETURN 
END 


sea  nest 

S3 


SUBROUTINE  INSERT! 2,  U.  I , NF.  ML ) 
SUBROUTINE  INSERT!  Z.  U.  I.  NF,  Nl  ) 


24 

S3 

26 

27 

20 

29 

30 

31 

32 

33 

34 
33 

36 

37 

38 

39 

40 

41 

42 

43 

44 
43 

46 

47 

48 

49 

30 

31 

32 

33 


30 

C 

C 

C 

C 


SUBROUTINE  USED  BY  SUBROUTINE  INSORT 
INSERT  1  NO  U< I ) 

(Z.W)ARE  DATA  PAIRS 

NR  IS  FIRST  ELEMENT  OF  ARRAY  ELIOIBLE  FOR  INSERTION 
NL  IS  LAST  ELEMENT  OF  ARRAY  ELIOIBLE  FOR  INSERTION 


DIMENSION  Z( 1). W< 1 > 


USINO  COMPILER  TIME  DOMAIN 
COMPUIATION  DONE  IN  PLACE 


J«I-l 
Y«W! I ) 

x-zd) 

IF t  Y  .  9E. 


W<U)  )  00  TO  100 


MOVES  U(J>  ONE  SPACE  US  AS 
THE  LEFT  OF  W!J> 


UK  Jfl)-U<  J) 
Z<J*1)-Z< J) 
ISTOP-NF-1 
J-J-t 


Y  WILL  BE  INSERTED  TO 


100 


IF  ( J  EQ. 
00  TO  30 
W!  J«-l  )-Y 
Z! J*1 >-X 
RETURN 
END 


I STOP)  CO  TO  100 


IFTRAN  STATISTICS 
33  CARDS  READ 
O  ERROR (S)  FOUND 
SEO  NEST 


SUBROUTINE  MXINVRCA;  B-  IBAD) 

SUBROUTINE  MXINVRIA,  B.  IBAD) 

C  SUBROUTINE  MXINVR  FACTORS  MATRIX  A  INTO  LU  CROUT  DECOMPOS I T I TON 
C  AND  COMPUTES  ITS  INVERSE;  MATRIX  B. 

IMPLICIT  DOUBLE  PRECISION  <A-H.O-Z> 


6 

7 

a 

9 

10 
1  I 
12 

13 

14 

13 

16 

17 

18 
19 


DIMENSION  A<3.  3).  B<3.  3) 


FOR ( I “ l  TO  3) 

FOR  <  J*1  TO  3) 

D(I.J)«0, 

IF (  I  .  EO.  J)  B  C  I.  J>< 
END  FOR 
END  FOR 


I8AD-0 

FORIK-l  TO  3) 

FOR <  I  "K  TO  3! 

DINN1-0  DO 
F0R1L-1  TO  ( K— 1 ) ) 


20 

3 

DINN1*DINN1+A< I, L)*A(L.  K> 

21 

2 

END  FOR 

22 

2 

A<  I, K)*A< I, K1-DINN1 

23 

1 

END  FOR 

24 

1 

C 

*  TEST  FOR  SINGULARITY. 

25 

1 

IF(  A  ( K<  K  > .  NE.  0  > 

26 

2 

FOR ( J*(K+1 >  TO  3) 

27 

3 

DINN2=0. DO 

28 

3 

FOR (L=l  TO  < K— 1 >  > 

29 

4 

DINN2»DINN2+A(K,  L>*A<L» U> 

30 

3 

END  FOR 

31 

3 

A(K, J)-DINN2)/A(K,  K) 

32 

2 

END  FOR 

33 

1 

ELSE 

34 

2 

IBAD»1 

35 

2 

K=4 

36 

1 

END  IF 

37 

END  FOR 

38 

c 

39 

IF(  IDAD.  NE.  1  ) 

40 

1 

FOR  ( M=1  TO  3) 

41 

2 

FOR (1=1  TO  3) 

42 

3 

D INN1=0.  DO 

43 

3 

FOR (K=l  TO  < I  —  1 >  > 

44 

4 

DINN1=DINN1+A< I,  K>*B(K,  M> 

>EQ 

NEST 

SUBROUTINE  MX INVR < A, B .  I B AD > 

45 

3 

END  FOR 

46 

3 

B< I, M)=<B< I,  MJ-DINN1 >/A< I,  I ) 

47 

2 

END  FOR 

48 

2 

FORC I»3  TO  1  STEP  -1) 

47 

3 

D INN 1=0.  DO 

50 

3 

FOR  ( J-  1  +  1  )  TO  3) 

51 

4 

DINN1=DINNH-A(  I,  U  )  *B  (  J,  M) 

52 

3 

END  FOR 

53 

3 

B(  I, M)=B( I, M1-DINN1 

54 

2 

END  FOR 

35 

1 

END  FOR 

56 

ELSE 

37 

1 

RETURN 

58 

END  IF 

59 

c 

60 

RETURN 

61 

END 

IFTRAN  STATISTICS 
61  CARDS  READ 
0  ERROR (S)  FOUND 


f  - —  -  ■  ““ — — “ - 

— _■ — - _ 

- 

- - 

Sj 

4 

<0001 ) 

SUBROUTINE  PDFIT  <N,  SDBS,  PS,  THS  ) 

<0002) 

C 

.  j 

<0003) 

c 

<0004) 

c 

jj 

<0005) 

c 

J 

<0006) 

IMPLICIT  DOUBLE  PRECISION  <A-H,  0-Z > 

'  -7‘|| 

<0007) 

c 

<0008) 

INTEGER  FILE6<3) , CODE 

J 

<0009) 

REAL  CASES,  PFAS.  THS (3) , SNR3( 12),  PCS< 12) ,  PDFS< 12) , 

(0010) 

1  ERRORS  ( 12) ,  ERRMS,  SDERRS,  SEC1,  SEC2-  SECS-  SES,  SDBS<  12),  C’S(  12) 

J 

<001  1  ) 

c 

<0012) 

DIMENSION  PD< 12) - ED( 12) , WK<3),  I STOP < 12),  ST< 12) , ERROR < 12) 

_r 

<0013) 

DIMENSION  SNR  < 12 > , PDF < 1 2 ) , TH < 3 ) ,  DEL < 3 )  -  DINC (3) ,  W<3,  12) 

<0014) 

DIMENSION  IP  < 3 )  -  STP  < 3 ) 

v. 

<0015) 

c 

<0016) 

COMMON/COM/SNR,  PDF,  TH,  DEL,  DINC,  W,  M,  KS,  JV 

i 

<0017) 

c 

i 

<0018) 

DATA  WK/O.  IDO,  0.  IDO,  0.  IDO/ 

<0019) 

DATA  ST/O.  00 1  DO,  0.  001D0,  0.  001D0,  0.  OOIDO,  0.  001D0,  0.  001D0, 

<0020) 

10.  OOIDO,  0.  OOIDO,  0.  OOIDO,  0.  OOIDO,  0.  OOIDO,  0.  OOIDO/ 

<0021  ) 

c 

~-4 

<0022) 

c 

INITIALIZATION. 

<0023) 

c 

<0024) 

DO  10  1  =  1, N 

<0025) 

PD  < I ) =DBLE  <  PS  < I ) ) 

<0026) 

SNR  < I >  =DBLE  < 10.  **  <  SDBS  < I )*.  1 ) ) 

<0027) 

10 

CONTINUE 

■j 

(0028) 

TH< 1 )=1.  D0/D3LE<  THS< 1 ) > 

! 

(0029) 

TH ( 2 ) =DBLE ( THS ( 2 ) ) 

(0030) 

TH  <  3 ) =D3LE ( THS  (  3 ) ) 

(0031 ) 

L=5 

<0032) 

M=0 

(0033) 

IP<1 )=0 

<0034) 

KS=0 

t 

(0033) 

IP  <  2 ) =0 

<0036) 

IPO)  =0 

<0037) 

ISS=0 

<0038) 

DEL (  1  )  =0.  000 IDO 

(0039) 

DEL  (  2  )  =0.  000 IDO 

J 

<0040) 

DEL  (  3  )  =0.  000 1  DO 

<0041 ) 

STP ( 1 )=0.  000001 DO 

'0042) 

STP ( 2 ) =0.  000001 DO 

(0043) 

STP  (  3 )  =0.  00000 IDO 

'1 

(0044) 

ISUM=0 

(0045) 

DMIN=1. OD-2 

(0046) 

WEIGH=1.  ODO 

(0047) 

SK= 1 .  ODO 

(0048) 

DO  20  1=1, 12 

i 

(0049) 

I  STOP  (  I  )  =0 

<0050) 

20 

CONTINUE 

<0051  ) 

JV=N 

<0 052) 

c 

(0050) 

c 

COMPUTE  INITIAL  DECUPLING  MATRIX. 

(0054) 

CALL  V3AR 

(0055) 

c 

( 0056 ) 

CALL  DCPLE 

"i 

(0057) 

KSMAX=500 

(0058) 

DO  160  KS=1, KSMAX 

(0059) 

M=M+1 

(0060) 

CALL  VBAR 

<0061 ) 

c 

(0062) 

c 

COMPUTE  DECUPLING  MATRIX 

(0063) 

c 

7 

(0064) 

IF ( M.  GT.  L)  CALL  DCPLE 

3 

<0065) 

c 

(0066) 

c 

FORM  ERRORS. 

(0067) 

c 

(0068) 

DO  110  1  =  1,  JV 

-d 

(0069) 

ED ( I )  =  ( PDF  < I ) -PD ( I ) >*WEIGH 

(0070) 

IF  (  DABS  (  ED  (  I  )  )  .  LE,  ST  (  I  )  )ISTOP(  I  >  =  1 

(0071 ) 

ISUM=  I SUM+ I  STOP ( I > 

(0072) 

110  CONTINUE 

= 

(0073) 

(0074) 

(0075) 

(0076) 

(0077) 

(0078) 

(0079) 

(ooao) 
(ooai ) 
(0082) 
(0083) 
(0084) 
(0085) 
(0086) 
(0087) 
( 0088 ) 
(0089) 
(0090) 
(0091  ) 
(0092) 
(0093) 
(0094) 
(0095) 
(0096) 
(0097) 
( C098 ) 
(0099) 
(0100) 
(0101  ) 
(0102) 
(0103) 
(0104) 
(0105) 
(0106) 
(0107) 
(0108) 
(0109) 
(0110) 
(0111) 
(0112) 
(0113) 
(0114) 
(0115) 
(0116) 
(0117) 
(0118) 
(0119) 
(0120) 
(0121  ) 
(0122) 
(0123) 
(0124) 
(0125 ) 
(0126) 
(0127) 
(0128) 
(0129) 
(0130) 
(0131  ) 
(0132) 
(0133) 
(0134) 
(0135) 
(0136) 
(0137) 
(0138) 
(0139) 
(0140) 
(0141  ) 
(0142) 
(0143) 
(0144 ) 


C 

C  COMPUTE  WEIGHTED  ERRORS. 

C 

DO  130  1=1, 3 
DINC < I  )  =0.  ODO 
DO  120  K= 1 , JV 

DINC  < I >=DINC( I >+W( I,  K)*ED(K> 

120  CONTINUE 
C 

C  COMPUTE  UPDATED  PARAMETER  ESTIMATE. 

C 

TH( I ) =TH ( I ) -WK( I ) *DINC ( I ) /SK 
SS=WK( I >*DINC< I ) 

I F ( DABS ( SS  ) .  LE.  STP  ( I  )  )  IP  (  I  )»1 
ISS=ISS+IP( I ) 

C 

C  CHECK  IF  ESTIMATE  HAS  CONVERGED.  IF  YES,  SET  ISTOP(I) 

C  EQUAL  TO  1. 

c 

130  CONTINUE 

c 

C  HAVE  ALL  THREE  ESTIMATES  CONVERGED? 

C 

IF ( ISUM.  EQ.  JV)  GO  TO  170 
I F ( ISS.  EQ.  3)  GO  TO  170 
DO  140  1  =  1 ,  JV 
I  STOP ( I ) =0 
140  CONTINUE 

DO  150  1=1, 3 
IP ( I ) =0 

150  CONTINUE 
ISS=0 
I SUM=0 

160  CONTINUE 
170  CONTINUE 

I F ( ISUM  .  EQ.  JV)WRITE( 1,  1001 ) 

IF  ( ISS  .EQ.  3JWRITEU,  1002) 

I F ( KS  .  GE.  KSMAX ;WRITE( 1, 1003) 

S2=0 .  ODO 

SS=0.  ODO 

DO  ISO  1  =  1,  JV 

ERROR (  I  >*PD< I )-PDF(  I  > 

S2=S2+ERR0R i I ) **2 
SS=SS+ERRQR ( I ) 

180  CONTINUE 

ERRM=SS/DBLE ( FLOAT ( JV ) ) 

SDERR=DSQRT ( S2/DELE ( FLOAT ( JV ) ) -ERRM*ERRM ) 

C 

C  DOUBLE  PRECISION  TO  SINGLE  PRECISION  CONVERSION. 

C 

CASES=SNGL ( CASE ) 

PFAS=SNGL(PFA) 

DO  190  1  =  1, 3 
THSi I ) =SNGL  <  TH ( I )  ) 

190  CONTINUE 

THS (  1  >»1.  /THS ( 1  ) 

DO  200  1  =  1 ,  JV 

SNRS (  I  )  =SNGL ( SNR ( I ) ) 

PDS( I >=SNGL(PD( I ) ) 

PDFS ( I ) =SNGL ( PDF ( I ) > 

ERRORS ( I ) =3NGL ( ERROR ( I ) ) 

200  CONTINUE 

ERRMS=SNOL ( ERRM  > 

SDERRS=SNGL ( SDERR > 

C 

C  WRITE  TO  OUTPUT  FILE. 

C 

CALL  PDFOUT ( N> THS, SNRS-  PDS, PDFS, ERRORS, ERRMS, SDERRS) 
C 

RETURN 

1001  FORMAT (  'PDF  I T  TERMINATION  1,  SMALL  ERRORS') 

1002  FORMAT (  'PDF  I T  TERMINATION  2,  SMALL  COF.  CHANGE') 


(0145)  1003  F0RMAT(  'PDFIT:  TERMINATION  3.  MAX.  ITERATIONS') 

(0146)  END 


CASE 

CASES 

CODE 

DABS 

DBLE 

DCPLE 

DEL 

DINC 

DM  IN 

DSQRT 

ED 

ERRM 

ER  RMS 

ERROR 

ERRORS 

FILES 

FLOAT 


IF 
ISS 
IS'  IP 
i  sun 

JV 

K 

KS 

KSMAX 

L 

M 

N 

PD 

PDF 

PDFIT 

PDFOUT 

PDFS 

PDS 

PFA 

PFAS 

PS 

S3 

SDBS 

SDERR 

SDERR5 

SEC1 

SEC2 

SECS 

SES 

SK 

SNGL 

SNR 

SNR  3 

SS 

ST 

5  TP 

TH 

THS 


D  EXTERNAL 
D  EXTERNAL 
D  EXTERNAL 
D  /COM/ 

D  /COM/ 

D 

D  EXTERNAL 

D 

D 

R 

D 

R 

I 

D  EXTERNAL 


001600 

001604 

000000 

000000 

000000 

ocoooo 

000154 

000170 

001616 

000000 

000104 

001622 

001626 

000164 

000244 

000274 

000000 

001630 


I  /COM/ 


I  /COM/ 

I 

I 

I  /COM/ 

I  ARGUMENT 
D 

D  /COM/ 

D 

D  EXTERNAL 

R 

R 

D 

R 

R  ARGUMENT 
D 

R  ARGUMENT 

D 

R 

R 

R 

R 

R 

r> 

R  EXTERNAL 

D  /COM/ 

R 

D 

C 

D 

D  /COM/ 


000277 

001633 

000302 

001634 

000426 

001635 
000425 
001636 
001637 
0C0424 
000003 
0003 16 
000060 
OCOOOO 
000000 
000376 
000426 
001640 
001644 
000005 
001650 
000004 
001654 
001660 
000000 
000000 
000000 
000000 
n01662 
UOOOOO 

000000 
000456 
001666 
00'  >24 
O0O5O6 
000140 


VBAR 

W 

WEIGH 

WK 


R  ARGUMENT  000006 

D  EXTERNAL  OCOOOO 
D  /COM/  000204 
g  001672 

D  000010 

000565 

001466 

001512 

001541 

000764 


0123A 
0009S 
0008S 
0070 
0025 
0056 
00133 
0013S 
0045M 
0119 
0012S 
01  ISM 
0009S 
0012S 
0009S 
0008S 
0118 
0024M 
0070 
0086 
0114 
0131 
0014S 
0037M 
0012S 
0044M 
0016S 
0113 
0078M 
001 68 
0057M 
0031M 
0016S 
0001S 
0012S 
0013S 
00015 
0140 
0009S 
0009S 
0124A 
0009S 
00018 
Oil  1M 
000  IS 
01  19M 
0009S 
0009S 
0009S 
0009S 
0009S 
0047M 
0123 
0135 
0013S 
0009S 
008SM 
0012S 
0014S 
0013S 
000  IS 
0140A 
0054 
0013S 
0046M 
0012S 

0024 

0108 

0109 

0110 

0063 


0123M 


0064 


0069M 
0119 
0135M 
01 14M 
0133M 

0119 

0025 

0071 

0087 

0115 

0132 

0033M 

0087M 

0049M 

0071M 

0051M 

0118 

0079 

0034M 

0058 

0064 

0032M 

0024 

0025M 

0016S 


0132M 

0131M 

0124M 
0009S 
Oil  5M 
0009S 
0136A 
0136M 


0084 

0124 

0136 

0016S 

0130M 

0086A 

00191 

004  1 M 

0016S 

0009S 

0060 

0016S 

0069 

00181 

007  9 
01  *. 
014 
0143D 
007 2D 


0028 

0029 

0030 

0118 

0038M 

0039M 

0040M 

0077M 

0079M 

0084 

0085 

0070A 

0079 

0135A 

0140A 

0115 

0116 

0133A 

0140A 

0026 

0048M 

0049 

0068M 

0076M 

0077 

0079 

0084 

0098M 

0099 

0101M 

0102 

0116 

0133 

0125M 

0126 

0129M 

0030M 

0036M 

0086M 

0087 

0097 

0104M 

0109 

0070M 

0071 

0099M 

0096 

0105M 

0108 

0069 

0078 

0096 

0090 

0119 

0129 

OOSBM 

0110 

0110 

0059M 

0064 

0051 

0140A 

0069 

0114 

0131A 

0069 

0114 

0132A 

01  19 


0069 
0085 
01  13M 
0130 

0102M 


0108 


0140A 

0140A 


0025A 

0119 

0026 

0140A 


0126 

0026M 

0140A 

01  12M 

0070 

0042M 

0028M 

0028A 


0079 


0130 

0130A 


0131  0132 


01 16M  0113 


0043M 

0029M 

0029A 


0086 

0030M  0084M 

0030A  0126M 


0133 


0126A 

0128M 


0084  0085 
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140 
1  50 
1  AO 
1170 
”1  SO 
1 90 
20 
200 


0000  ERRORS 


001032 
001103 
001 126 
001141 
001152 
001160 
001244 
001346 
000700 
001430 


0078 

0076 

0098 

0101 

0058 

0096 

0113 

0125 

0040 

0129 
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0080D 

0092D 

0100D 

0103D 

01060 

0097  0107D 

01170 

0127D 

00500 

01340 


(0001  ) 
(0002) 
(0003) 
(0003) 
(0004) 
(0004) 
(0005) 
(0005) 
(0007) 
(0008) 
(0009) 
(0010) 
(OOll ) 
(0012) 
(0013) 
(0014) 
(0013) 
(0015) 
(0017) 
(0018) 
(0019) 
(0020) 
(0021 ) 
(0022) 
(0023) 
(0024) 
(0025) 
(0026) 
(0027) 
(0028) 
(0029) 
(0030) 
(0031 ) 
(0032) 
(0033) 
(0034) 
(0035) 
(0036) 
(0037) 
(0038) 
(0039) 
(0040) 
(0041 ) 
(0042) 
(0043) 
(0044) 
(0045) 
(0046) 
(0047) 
(0048) 
(0049) 
( 0030 i 
(0051 ) 
(0052) 
(0053) 
(0054 ) 
(0055) 
(0056) 
(0057) 
(0058) 
(0059) 
(0050) 
(0061 ) 
(0062) 
(0063) 
(0064 ) 
(0063) 
(0066) 
( 0067 ) 
(0068) 
(00o9) 
<00'7C) 


SUBROUTINE  PDFOUT ( N, TH,  SNR,  PD,  PDF,  ERROR,  ERRM, SDERR ) 
:OM>KEYS.  F  MNEMONIC  KEYS  FOR  FILE  SYSTEM  ( FTN ) 


C  SYSCOM>KEYS.  F  MNEMONIC  KEYS  FOR  FILE  SYSTEM 

NOLIST 

C  SYSCOM>ERRD.  F  MNEMONIC  CODES  FOR  FILE  SYSTEM  (T TN> 

NOLIST 

INTEGER  FILE6(6),  CODE 

DIMENSION  TH  (  3  )  #  SNR  U  > .  PDd  ) ,  PDF<  1  ) .  ERROR  <  1  > 


31  MAY,  197 
6  SEPT,  19 


ENTER  OUTPUT  FILE  NAME. 

CONTINUE 
WRITEd,  1200) 

READ  <  1,  1300)FIL£6 

SEARCH  TO  SEE  IF  FILE  ALREADY  EXISTS. 

CALL  SRCH®4 (K4EXST,  FILE6,  12,  2,  0,  ICODE) 

IF(  I  CODE  .  EQ.  E*FNTF)GO  TO  20 
WRITEd,  1006) 

READ ( 1 ,  *) I OK 
IF ( I OK  EQ.  1  >G0  TO  20 
GO  TO  10 
CONTINUE 

OPEN  OUTPUT  FILE. 

CALL  SRCH$i(K<WR IT ,  FILE6,  12,  2,  0, CODE) 

IF(  CODE.  NE.  0  )  CALL  ERRPR*  (K*NRTN,  CODE,  O,  O,  'PFDOUT  ' ,  6  ) 


C 

6666 


WRITEd,  6666) 

>6  FORMAT(  'ENTER:  CASE  PFA  NP  NB  '  /  ) 

READ  (  1,  *  )  CASE,  PFA,  NP,  NB 

GO  TO  (210,220,230,240),  NB 
)  WRITE<6,  1431 ) 

GO  TO  250 
)  WR  I TE  ( 6,  1432) 

GO  TO  250 
)  WRITE  (6,  1433) 

GO  TO  250 
)  WR  I  TE  <  6,  1434) 

)  WRITE  (6,  1440) 

WRITE16,  1460) 

WRITE  (6,  1500) 

WR I ^E  <6,  1600) 

WR  I  TE  ( 6,  1700)CASE,  PFA,  NP 
WR  I  TE  ( 6,  1800) 

WR  I  TE  (6,  1900) 

WR  I  TE  (6,  2000) 

WR  I  TE  ( 6,  2100) 

WR  I  TE  ( 6,  2200  )  TH  d  > ,  TH  (  2  )  ,  TH  (  3  ) 

WRITE (6,  2300) 

WR  I  TE  (6,  2400) 

WR  I  TE  (6,  2500) 

WR  I  TE  (  6,  2600) 

WRITE  (6,  2700)  <  I,  SNR  (  I  ) ,  PD  (  I  )  -  PDF  (  I  ) ,  ERROR  (  I  > ,  1*1,  N> 
WR  I  TE  ( 6,  2S00 )  ERRM,  SDERR 

REWIND  OUTPUT  FILE. 

REWIND  6 

CLOSE  OUTPUT  FILE. 

CALL  SRCH**(K*CLOS,  0,  0,  2,  0,  CODE) 

IF( CODE  EQ  0)  GO  TO  300 

CALL  ERRPR*(K*NRTN,  CODE,  0,  0,  'FDFOUT',6) 

}  CONTINUE 


(0071) 
(0073) 
(0073) 
(0074) 
(0073) 
(0076) 
(0077) 
(0078) 
(0079) 
( 0000 ) 
( 0001 ) 
(0002) 
(0083) 
(0004) 
(0003) 
<  0006 ) 
(00B7) 
(0088) 
(0089) 
(0090) 
(0091 ) 
(0092) 
(0093) 
(0094) 
(0093) 
1 0096 ) 
(0097) 
( 0098 ) 
( 0099 ) 
(0100) 


RETURN 

FORMATS, 

S§  "°t'  ”,r”  eLSt 

31  FORMAT  (  /  // /36X>  'RAVLEIOH  '  > 

2?  FORMAT! ////36X,  'LOO  -  NORMAL  -303') 

n  :fa?s BGLr-^3-,6DB,> 

l°o 

30  iP0^^1?JteI^t^  CAaE  '•  10x'  'PFA''  l2x- 

30  FORMAT  ( 13X<  ' - —  — — '.10X,  ',12X  ' _ 

DO  FORMAT(//22X,F3.  1,  UX,£11  1.  17X,  13)  1 

DO  rORMftT(31)(/:^4X'  ,PARAHETER  ESTATES') 

DO  FORMAT  (iex??---~~5IE? i'iJf  <1^*S22S7!  .2[q^  'POSTER  3 ' ) 

58  S8ja?{?«55i5,1W?*M4-7'w*M4-7‘ 

DO  FORMAT  (  42X.  - - - 

30  'S^R .'^i0X'  'PD''  lOX'  ,pDF'.  10X,  'ERROR') 

JO  FORMAT!  1BX>  ,  10X,  ' - lOX,  10X,  ' ,  lox.  — M 

DO  F0RMAT(17X,  I2.6X.E11,  4.  3X.F0.  5,4X,  FB  3.  6X.F8  3) 

DO  FORMAT!  ///19X  'MEAw  ERROR  •  ' ,  F7.  3,  12X,  'RMS  ERROR  -  ' ,  F7  3) 


rnnl  ?  001337 

CODE  I  001361 

|*bcqd  1  parameter 

E*BDAM  I  PARAMETER 
f*9£SV  I  PARAMETER 
EfBFTS  I  PARAMETER 
E*BKEY  I  PARAMETER 
E*BNAM  I  PARAMETER 
E*B0F  I  PARAMETER 
|*|PAR  I  PARAMETER 
E*BPAS  I  PARAMETER 
E*BSGN  I  PARAMETER 
|*3SUN  I  PARAMETER 
E*BTRN  I  PARAMETER 
E*BUFD  I  PARAMETER 
E*BUNT  I  PARAMETER 
E*DIRE  I  PARAMETER 
E*DISK  I  PARAMETER 
E*DKFL  I  PARAMETER 
E5DNS  I  PARAMETER 
ESONTE  I  PARAMETER 
E*DVIU  I  PARAMETER 
E*EOF  1  PARAMETEP 
E*EXST  I  PARAMETER 
ffPABT  I  PARAMETER 
|*PB3T  I  PARAMETER 
|*FDEL  I  PARAMETER 
E*FDFL  I  PARAMETER 
E*FIFC  I  PARAMETER 
E*FITB  I  PARAMETER 
E*FIUS  I  PARAMETER 

|*£NIF  1  parameter 

E*FNTS  I  PARAMETER 
E*FONC  I  PARAMETER 
|*^ru  I  PARAMETER 
E*IREM  I  PARAMETER 
IJNASS  I  PARAMETER 
EtNATT  I  PARAMETER 
|*NMLG  I  PARAMETER 
|*NPHA  1  PARAMETER 
EtNRIT  I  PARAMETER 
E3SNTIM  1  PARAMETER 
E$NTSD  I  PARAMETER 
E*NTUD  I  PARAMETER 
F.SNULL  I  PARAMETER 
EJOLDP  I  PARAMETER 
EiPTRM  I  PARAMETER 
£*PLON  I  PARAMETER 
ESRCOM  I  PARAMETER 
E*SDER  I  PARAMETER 
EtSEMO  I  PARAMETER 
ESSHUT  I  PARAMETER 
i*SUNO  I  PARAMETER 
E*TMRU  I  PARAMETER 
E*TMUL  I  PARAMETER 
ESUIUS  I  PARAMETER 
E*UNOP  I  PARAMETER 
ERRM  R  ARGUMENT  000011 
15555  R  ARGUMENT  000010 
ERRPR*  R  EXTERNAL  000000 
FILE6  I  000014 

Jrnnc  J  001362 

JCODE  I  001363 

*0K  1  001364 

K*ALLD  I  PARAMETER 
K*CACC  I  PARAMETER 
K*CLOS  I  PARAMETER 
K*CGNV  I  PARAMETER 
K*CURR  I  PARAMETER 
KiDELE  I  PARAMETER 
K*DMPD  I  PARAMETER 
K4DTIM  I  PARAMETER 


0033M  0046 

OOoll  0026A  O027A  0067A  0068 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0C04S 

0004S 

0004S 

0004S  0017 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

0004S 

000 IS  0059 
0001S  0006S  003fl 

0027  0069 

0038M  °012f1  0016A  0026A 

0016A  0017 

0019M  0020 

0003S 

0003S 

0003S  0067 

0003S 

0003S 

0003S 

0003S 

0003S 


006VA 


1 


K*ENTR 

I 

000000 

0003S 

KSEXST 

I 

PARAMETER 

0003S 

0016 

K*FREE 

I 

PARAMETER 

0003S 

K*FULL 

I 

PARAMETER 

0003S 

KJGOND 

I 

PARAMETER 

0003S 

K*OPOS 

1 

PARAMETER 

0003S 

K*HOKE 

I 

PARAMETER 

0003S 

KJICUR 

I 

PARAMETER 

0003S 

KSIMFD 

I 

PARAMETER 

0003S 

►<*IRTN 

I 

PARAMETER 

0003S 

K*1SEG 

I 

PARAMETER 

0003S 

K*IUFD 

I 

PARAMETER 

0003S 

K*MENT 

1 

000000 

0003S 

K*MS1  7. 

I 

PARAMETER 

0003S 

K*MVNT 

I 

PARAMETER 

0003S 

K*NAME 

I 

PARAMETER 

0003S 

K*NDAM 

I 

PARAMETER 

0003S 

KiNRTN 

I 

PARAMETER 

0003S 

0027 

0069 

K*NSAM 

I 

PARAMETER 

0003S 

K4NSCD 

I 

PARAMETER 

0003S 

K*NSGS 

I 

PARAMETER 

0003S 

K*POSA 

I 

PARAMETER 

0003S 

K*POSN 

I 

PARAMETER 

0003S 

K*POSR 

I 

PARAMETER 

0003S 

K*PREA 

I 

PARAMETER 

0003S 

K$PRER 

I 

PARAMETER 

0003S 

K*PROT 

I 

PARAMETER 

0003S 

K$RDWR 

I 

PARAMETER 

0003S 

K*READ 

I 

PARAMETER 

0003S 

K*RPOS 

I 

PARAMETER 

0003S 

K*RS'JB 

I 

PARAMETER 

0003S 

KiRWLK 

1 

PARAMETER 

0003S 

K*SENT 

I 

000000 

0003S 

K*SETC 

I 

PARAMETER 

0003S 

K*SETH 

I 

PARAMETER 

0003S 

K*SPOS 

I 

PARAMETER 

0003S 

K*SRTN 

I 

PARAMETER 

0003S 

K*TRNC 

I 

PARAMETER 

0003S 

K*UPOS 

I 

PARAMETER 

0003S 

K*URIT 

I 

PARAMETER 

0003S 

0026 

N 

I 

ARGUMENT  000003 

000  IS 

0058 

NB 

I 

001365 

0033M 

0035 

NP 

I 

001366 

0033M 

0048 

PD 

R 

ARGUMENT  000006 

0001S 

00065 

0058 

PDF 

R 

ARGUMENT  000007 

000  IS 

0006S 

0058 

PDFOUT 

R 

OCOOOO 

0001S 

PFA 

R 

001367 

0033M 

0048 

SDERR 

R 

ARGUMENT  000012 

0001S 

0059 

SNR 

R 

ARGUMENT  000005 

000  IS 

0006S 

0058 

SRCH** 

R 

EXTERNAL  000000 

0016 

0026 

0067 

TH 

R 

ARGUMENT  000004 

000  IS 

00068 

0053 

10 

000022 

0010D 

0021 

1006 

000457 

0018 

0075D 

”1200 

000515 

0011 

0076D 

“1300 

000537 

0012 

0077D 

“1431 

000543 

0036 

007SD 

1432 

000556 

0038 

0079D 

“1433 

000576 

0040 

0080D 

“1434 

000616 

0042 

0081  D 

“1440 

000634 

0044 

0082D 

“1460 

000653 

0045 

0083D 

“1  500 

000670 

0046 

0084D 

“1600 

000730 

0047 

0086D 

“1700 

000767 

0048 

0087D 

1800 

001007 

0049 

0088D 

“1  900 

001030 

0050 

00S9D 

“20 

000071 

0017 

0020 

0022D 

“2000 

001047 

0051 

0090D 

“210 

000165 

0035 

0036D 

“2100 

001 105 

0052 

009  ID 

“220 

000172 

0035 

0038D 

2200 

001142 

0053 

0092D 

”230 

000177 

0033 

0040D 

"2300 

001163 

0034 

0093D 

"240 

000204 

0033 

0042D 

"2400 

001174 

0053 

0094D 

"250 

000210 

0037 

0039 

"2500 

001203 

0036 

009SD 

2600 

001236 

0037 

0096D 

"2700 

001270 

0038 

0097D 

"2800 

001316 

0039 

0098D 

"300 

000453 

0068 

0070D 

"6666 

000120 

0030 

0032D 
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0044D 


(0001  ) 
(0002) 
(0003) 
(0004) 
(0005) 
(0006) 
(0007) 

( 0003 ) 
(0009) 
(0010) 
(OOl  i  ) 
(0012) 
(0013) 
(0014) 
(0015) 
(0016) 
(0017) 
(0018) 
(0019) 
( 0020 ) 
(0021  ) 
(0022) 
(0023) 
(0024) 
(0025) 
( 0026 ) 
(0027) 
( 0028 ) 
( 0029 ) 
( 0030 ) 
(0031  ) 
(0032) 
(0033) 
(0034) 
(0035) 
(0036) 
(0037) 
( 0038 ) 
(0039) 
(0040) 
(0041  ) 
(0042) 
(0043) 
(0044) 
(0045) 
(0046) 
(0047) 
(0048) 
(0049) 
(0050) 
(0051 ) 
(0052) 
(0053) 
(0054) 
(0055) 
(0056) 
(0057) 
(0058) 
( 005° ' 
( 0060 ) 
(0061 ) 
(0062) 
(0063) 
(0064) 
(0065) 
(0066) 
(0067) 
(0068 ) 
(0069) 
(0070) 
(0071  ) 
(0072) 

(0073) 


C 

C 

c 


10 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


SUBROUTINE  PTDATA ( N. X,  Y ) 

DIMENSION  X( 1 >. Y( 1 >. LABLEX( 16). LABLEY < 16) 

COMMON  /PTCOM/XORIG, YORIG, XMAX. XMIN, YMIN. SCALEX, SCALEYi 
*S I ZEX 

FIND  MAX  S<  MIN  VALUES  OF  X  5c  Y. 

XMIN-1.  E32 
YMIN-1.  E32 
XMAX— XMIN 
YMAX— YMIN 


DO  10  I-l.N 

IF < X(  I  ) 

.  GT. 

XMAX)XMAX-X(I) 

IF( X ( I  > 

.  LT. 

XMIN) XMIN=X ( I ) 

IF <  Y( I ) 

.  GT. 

YMAX ) YMAX=Y( I ) 

IF<  Y( I > 

.  LT. 

YMIN) YMIN-YC I ) 

100 

C 

C 

C 

c 

1001 

1002 

1003 

1004 
C 


CONTINUE 

INPUT  PLOT  SIZE. 

WRITE< 1. 1001 > 

READ ( 1 1  *>SIZEX.  SIZEY 

RE-ADJUST  XMAX-XMIN,  YMAX-YMIN  TO  MAKE  A  NICE  PLOT 

CALL  SCALF( XMAX.  XMIN, DELX, ANTICX ) 

CALL  SCALE (YMAX. YMIN, DELY,  ANTICY) 

COMPUTE  SCALE  FACTORS 

SC ALEX = ( XMAX— XMIN ) /SI ZEX 
SCALEY=( YMAX-YMIN) /SIZEY 
H=. 022*SIZEX 

SET  ORIGIN 

X0RIG=2.  0 
YORIG-1.  5 

CALL  PLOT(XORIG,  YORIG,  -3) 

DRAW  GRID  5<  LABEL. 

NT I CX* IF IX (ANTICX ) 

NTICY-IFIX(ANTICY) 

CALL  BORDER  (SI  ZEX,  SIZEY,  NTICX.  NTICY,  1,  1,  H) 

WRITE(1, 1002) 

READ  < 1 1  1003 ) LABLEX 
WRITE( 1, 1004) 

READ ( 1 , 1003 ) LABLEY 

CALL  AXLABL ( SI ZEX . SIZEY, H,  XMIN,  DELX,  XMAX,  1,  1, YMIN, DELY, 
*YMAX ,  1,1, LABLEX, 32, LABLEY, 32  > 

PLOT  DATA 

DO  100  1  =  1,  N 

XP»(X< I ) -XMIN) /SCALEX 
YP»( Y( I )-YMIN) /SCALEY 
CALL  PLOT  <  XP ,  YP, 3) 

H2=H 

CALL  CSYMBL<XP, YP, H2, 2,  O.  , -2) 

CONTINUE 


RE-SET  ORIGIN. 

CALL  PLOT(-XORIG, -YORIG, -3) 

FORMAT (  'ENTER  XDIM( IN) ,  YDIM( IN)  ' ) 

FORMAT(  'ENTER  X  AXIS  LABEL,  FORMAT  32A1  ' ) 


FORMAT ( 32A1 ) 

FORMAT (  'ENTER  YAXIS  LABEL, 

RETURN 


FORMAT  32A 1 ' ) 


END 


ANT ICX 

R 

000333 

0026A 

0043 

ANTICY 

R 

000535 

0027A 

0044 

AXLABl 

R 

EXTERNAL 

000000 

0050 

BORDER 

R 

EXTERNAL 

000000 

0043 

CSYMBL 

R 

EXTERNAL 

000000 

0060 

DEUX 

R 

000337 

0026A 

0030A 

DELY 

R 

000341 

0027A 

0030A 

H 

R 

000343 

0033M 

004  3A 

0030A 

0039 

H2 

R 

000343 

00S9M 

0060A 

I 

I 

000347 

0012M 

0057 

0013 

0014 

0015 

0016 

0033M 

0036 

IF  I  X 

I 

EXTERNAL 

000000 

0043 

0044 

LABLEX 

I 

000007 

00025 

0047M 

0030A 

LABLEY 

I 

000027 

0002S 

0049M 

0030A 

N 

I 

ARGUMENT 

000003 

000  IS 

0012 

0055 

NTICX 

I 

000352 

0043M 

0043A 

NTICY 

I 

000553 

0044M 

0043A 

PLOT 

R 

EXTERNAL 

000000 

0039 

0058 

0065 

PTDATA 

R 

OOCOOO 

000  IS 

SCALE 

R 

EXTERNAL 

000000 

0026 

0027 

SCALEX 

R 

/PTCOM/ 

000012 

0003S 

0031M 

0036 

SCALEY 

R 

/PTCCM/ 

000014 

0003S 

0032M 

0037 

SIZEX 

R 

/PTCOM/ 

000016 

0003S 

0022M 

0031 

0033 

0045A 

0050A 

SIZEY 

R 

000360 

0022M 

0032 

0043A 

0030A 

X 

R 

ARGUMENT 

000004 

0001S 

0002S 

0013 

0014 

0056 

XMAX 

R 

/PTCOM/ 

000004 

0003S 

0010M 

0013M 

0026A 

0031 

0050A 

XMIN 

R 

/PTCOM/ 

000006 

0003S 

0056 

0008M 

0010 

0014M 

0026A 

0031 

0030A 

XORIG 

R 

/PTCOM/ 

000000 

0003S 

0037M 

0039A 

0063 

XP 

R 

000562 

0036M 

0058A 

0060A 

Y 

R 

ARGUMENT 

000005 

0001S 

0002S 

0015 

0016 

0057 

YMAX 

R 

000564 

001 1M 

0015M 

0027A 

0032 

0050A 

YMIN 

R 

/PTCOM/ 

000010 

0003S 

0057 

0009M 

001 1 

0016M 

0027A 

0032 

0030A 

YORIG 

R 

/PTCOM/ 

000002 

0003S 

0038M 

0039A 

0065 

YP 

R 

000566 

0037M 

0058A 

0060A 

10 

000144 

0012 

0017D 

100 

000413 

0053 

006  ID 

lOOl 

000441 

0021 

0067D 

1002 

000460 

0046 

0068D 

“1003 

000503 

0047 

0049 

0069D 

"1004 

000507 

004S 

0070D 
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(0001 ) 
(0002) 
(0003) 
(0004) 
(0005) 
( 0006 ) 
(0007) 
(0008) 
( 0009  > 
(0010) 
(0011  ) 
(0012) 
(0013) 
(0014) 
(0015) 
(0016) 
(0017) 
(0018) 
(0019) 
( 0020 ) 
(0021  ) 
( 0022 ) 
(0023) 
(0024) 
(0023) 
( 0026 ) 
(0027) 
(0028) 
( 0029 ) 
(0030) 
(0031  ) 
( 0032 ) 
(0033) 
(0034) 


SUBROUTINE  PTQFIT(AK) 

C 

DIMENSION  AK<3> 

COMMON  /PTCOM/XORIG,  YORIG,  XMAX,  XMIN,  YMIN,  SCALEX,  SCAl.EY, 
*SIZEX 
C 

PF(A.  B.  C,  S>«1.  -  ( 1. +(A#S>**8 >**(-C> 

C 

C  SET  ORIGIN. 

C 

CALL  PLOTCXORIG,  YORIG, -3) 

C 

C  COMPUTE  VALUES  Sc  PLOT. 

C 

J*>1 

X»XMAX 

DELX=»SI  ZEX/100. 

10  IF ( X  . LT.  XMIN)GO  TO  100 

Y=PF <  AK  < 1 ) ,  AK ( 2  > ,  AK ( 3  > .  X  > 

XP*( X-XMIN) /SCALEX 
YP*( Y-YMIN) /SCALEY 
IF ( J  .  EQ.  1 ) CALL  PLOT( XP, YP,  3) 

CALL  PLOT ( XP , YP, 2) 

J*2 

X»X-DELX 
GO  TO  10 
100  CONTINUE 

C 

C  RESET  ORIGIN. 

C 

CALL  PLOT(-XORIG, -YORIG, -3) 

C 

RETURN 

END 


A 

R 

000000 

0007 

AK 

R 

ARGUMENT 

000003 

0001S 

0003S 

0019A 

B 

R 

000000 

0007 

C 

R 

000000 

0007 

DEL* 

R 

000201 

0017M 

0025 

J 

I 

000203 

0013M 

0022 

0024M 

PE 

R 

000005 

0007S 

0019 

PLOT 

R 

EXTERNAL 

000000 

0011 

0022 

0023 

PTQFIT 

R 

000000 

0001S 

S 

R 

000000 

0007 

SCALEX 

R 

/P TCOM/ 

000012 

0004S 

0020 

SCALEY 

R 

/P TCCM/ 

000014 

0004S 

0021 

SIZEX 

R 

/P TCOM/ 

000016 

0004S 

0017 

X 

R 

000212 

0016M 

0018 

0019A 

XMAX 

R 

/PTCOM/ 

000004 

0004S 

0016 

XMIN 

R 

/P TCOM/ 

000006 

0004S 

0018 

0020 

XORIG 

R 

/PTCOM/ 

000000 

0004S 

001 1A 

0031 

XP 

R 

000214 

0020M 

0022A 

0023A 

Y 

R 

0C0216 

0019M 

0021 

YMIN 

R 

/PTCOM/ 

000010 

0004S 

0021 

YORIG 

R 

/PTCOM/ 

000002 

0004S 

001 1A 

0031 

YP 

R 

000220 

002 1M 

0022A 

0023A 

10 

000062 

001BD 

0026 

_1 00 

000157 

0018 

0027D 
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0031 


0020 


002  ISM 


(0001  ) 

SUBROUTINE  SCALE (RMAX, RMIN.  RINC,  ANTIC  > 

(0002) 

RMIN=0. 

(0003) 

IF ( RMAX  .  GT.  1.  ) GO  TO  5 

(0004) 

RMAX»1. 

(0003) 

GO  TO  100 

( 0006 ) 

5 

IF ( RMAX  .  GT.  10.  >G0  TO  10 

(0007) 

RMAX-10. 

( 0008 ) 

GO  TO  100 

(0009) 

10 

IF(RMAX  .  GT.  20.  >G0  TO  20 

(0010) 

RMAX*20. 

(0011  ) 

GO  TO  100 

(0012) 

20 

IF ( RMAX  .  GT.  SO.  >G0  TO  30 

(0013) 

RMAX*30. 

(0014) 

GO  TO  100 

(0015) 

30 

IF (RMAX  .  GT  100.  >G0  TO  40 

(0016) 

RMAX=  1 00. 

(0017) 

GO  TO  100 

(0018) 

40 

IF (RMAX  .  GT  130.  >G0  TO  30 

(0019) 

RMAX  =  1  50. 

( 0020 ) 

GO  TO  100 

(0021  ) 

50 

I F ( RMAX  .  GT.  200.  >G0  TO  60 

(0022) 

RMAX=200. 

(0023) 

GO  TO  100 

(0024 ) 

60 

I F ( RMAX  GT  230.  )G0  TO  70 

(0025) 

RMAX*250. 

( 0026 ) 

GO  TO  100 

(0027) 

70 

IF ( RMAX  .  GT.  300.  >G0  TO  80 

(0028) 

RMAX=300. 

(0029) 

GO  TO  100 

(0030) 

80 

RMAX  =  500. 

(0031 ) 

100 

CONTINUE 

(0032) 

R INC=( RMAX— RMIN ) / 10. 

(0033) 

ANTIC=10. 

(0034) 

RETURN 

(0035) 

END 

n_ 

. 

a 

3 

* 

I 

?  I' 

ANTIC 

R 

ARGUMENT 

000006 

0001S 

0033M 

• 

'*  ?  ;  . 

RINC 

R 

ARGUMENT 

000005 

000  IS 

0032M 

4 

i 

-  -  - 

RMAX 

R 

ARGUMENT 

000003 

000  IS 

0003 

0004M 

0006 

0007M 

0009 

0010M 

- 

0012 

0013M 

0015 

0016M 

0018 

0019M 

0021 

: 

0022M 

0024 

0025M 

0027 

0028M 

0030M 

0032 

RMIN 

R 

ARGUMENT 

000004 

00U1S 

0002M 

0032 

SCALE 

R 

000000 

000  IS 

% 

E- 

10 

000037 

0006 

0009D 

• 

“l  00 

000151 

0005 

0008 

0011 

0014 

0017 

0020 

0023 

- 

0026 

0029 

0031D 

20 

000051 

0009 

001 2D 

- 

“30 

000063 

0012 

001  5D 

- 

“40 

000075 

0015 

0018D 

4 

“5 

000025 

0003 

0006D 

. 

“50 

000107 

0018 

002  ID 

“60 

000121 

0021 

0024D 

“70 

000133 

0024 

0027D 

■  * 

| 

“BO 

000145 

0027 

0030D 

U 

I 

; 
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SUBROUTINE  QFIN(N, SDB, P > 


(0001 ) 

(0002) 

C 

(0003) 

(0003) 

(0004) 

(0005) 

C  SY! 

(0006) 

C 

(0007) 

C 

(0008) 

(0009) 

(ooio; 

C 

(0011 ) 

C 

(0012) 

C 

(0013) 

(0014) 

(0015) 

C 

(0016) 

C 

(0017) 

C 

(001B) 

(0019) 

( 0020 ) 
(0021 ) 
(0022) 

C 

( 0023 ) 

50 

( 0024 ) 

C 

(0025) 

C 

(0026) 

(0027) 

(0028) 

C 

(0029) 

(0030) 

C 

(0031 ) 

C 

( 0032 ) 

1001 

(0033) 

1002 

(0034) 

1003 

(0035) 

(0036) 

c 

:OM>KEYS. F  MNEMONIC  KEYS  FOR  FILE  SYSTEM  ( FTN ) 

NOLIST 

DIMENSION  SDB ( 1 ) . P  < 1 » 

DIMENSION  NAMEF ( 6 ) 

READ  INPUT  FILE  NAME. 

WR I TE ( 1 , 1001 ) 

READ  < 1 ,  1002) NAMEF 

OPEN  FILE  FOR  READING. 

CALL  SRCH**(K*READ,  NAMEF,  12,  1,0,  ICODE) 

Ir  < I  CODE  .  NE.  0)CALL  ERRPR*  (  K*NRTN,  I  CODE,  0,  0,  O.  0  ) 

READ  SDB, P  VALUES. 

READ  <5,  *  >  N 

IF ( N  .  GT  50>WRITE< 1,  1003) 

DO  50  1*1,  N 

READO,  *  )  SDB  (  I  ) .  P  <  I  > 

CONTINUE 

CLOSE  INPUT  FILE. 

CALL  SRCH**<K*CLOS,  NAMEF,  12,  1,0,  ICODE) 

IF  (  ICODE  .  NE.  OCALL  ERRPR*  (  K*NRTN,  ICODE,  O,  O,  O,  O  > 


31  MAN 


RETURN 

FORMAT (  ENTER 
FORMAT (AA2) 
FORMAT (  '  INPUT: 

END 


INPUT  FILE  NAME,  FORMAT  6A2 ' ) 

****  ERROR,  TOO  MANY  DATA  PAIRS  FOR  STORAGE  ** 


i—===.^F=  -  -  *  ‘  -  .  ^ 

— - -  - 

'  i 

F 

I  •  <“ 

1 

| 

ERRPR* 

R 

EXTERNAL  000000 

0015 

0028 

- 

fi 

I 

I 

0C0244 

002 1M 

0022 

1 

ICODE 

I 

000245 

0014A 

0015A 

0027A 

002SA 

r 

KJALLD 

I 

PARAMETER 

0003S 

K4CACC 

I 

PARAMETER 

0003S 

1 

K4CL0S 

I 

PARAMETER 

0003S 

0027 

fr  r  ; 

KS CONV 

I 

PARAMETER 

0003S 

\ 

r  .  , 

K4CURR 

I 

PARAMETER 

0003S 

l 

K4DELE 

I 

PARAMETER 

0003S 

K4 DMPB 

I 

PARAMETER 

0003S 

t 

K4DTIM 

I 

PARAMETER 

0003S 

K4ENTR 

I 

000000 

0003S 

K4EXST 

I 

PARAMETER 

0003S 

L 

K4FREE 

I 

PARAMETER 

0003S 

f 

p. 

K4FULL 

I 

PARAMETER 

0003S 

! 

K4GOND 

I 

PARAMETER 

0003S 

7 

KiGPOS 

I 

PARAMETER 

0003S 

K4HOME 

I 

PARAMETER 

0003S 

K2ICUR 

I 

PARAMETER 

0003S 

K4 I MFD 

I 

PARAMETER 

0003S 

K$ IRTN 

I 

PARAMETER 

0003S 

Kt ISEG 

I 

PARAMETER 

0003S 

K$ IUFD 

I 

PARAMETER 

0003S 

K4MENT 

I 

000000 

0003S 

K4MSI7 

I 

PARAMETER 

0003S 

i 

K4MVNT 

I 

PARAMETER 

0003S 

K4NAME 

I 

PARAMETER 

0003S 

fe 

K4NDAM 

I 

PARAMETER 

0003S 

K4NRTN 

I 

PARAMETER 

0003S 

0015 

0028 

K4NSAM 

I 

PARAMETER 

0003S 

K4NSGD 

I 

PARAMETER 

0003S 

KtMSGS 

I 

PARAMETER 

0003S 

KSPOSA 

I 

PARAMETER 

0003S 

K4POSN 

1 

PARAMETER 

0003S 

K4P0SR 

I 

PARAMETER 

0003S 

- 

K4PREA 

I 

PARAMETER 

0003S 

K4PRER 

I 

PARAMETER 

0003S 

K4PR0T 

I 

PARAMETER 

0003S 

K4RDWR 

I 

PARAMETER 

0003S 

K4READ 

I 

PARAMETER 

0003S 

0014 

K4RPQS 

I 

PARAMETER 

0003S 

s. 

K4RSUB 

I 

PARAMETER 

0003S 

r 

K4RWLK 

I 

PARAMETER 

0003S 

1 

K4SENT 

I 

OCOOOO 

0003S 

1 

K4SETC 

I 

PARAMETER 

0003S 

I 

K4SETH 

I 

PARAMETER 

0003S 

I 

K4SPOS 

I 

PARAMETER 

0003S 

1 

K4SRTN 

I 

PARAMETER 

0003S 

f 

K4TRNC 

I 

PARAMETER 

0003S 

i 

f 

K4UPOS 

I 

PARAMETER 

0003S 

p 

K4WRIT 

I 

PARAMETER 

0003S 

N 

I 

ARGUMENT  000003 

000  IS 

0019M 

0020 

0021 

NAMEF 

I 

000007 

0005S 

0010M 

0014A 

0027A 

P 

R 

ARGUMENT  000005 

000  IS 

0004S 

0022M 

QF  IN 

R 

000000 

000  IS 

SDB 

R 

ARGUMENT  000004 

000  IS 

0004S 

0022M 

SRCH9* 

R 

EXTERNAL  OOOOOO 

0014 

0027 

1001 

0C0154 

0009 

0032D 

t 

1002 

000177 

0010 

0032D 

t 

1003 

000203 

0020 

0034D 

1 

50 

000122 

0021 

00230 

r 

t 

i 

1 

b 

F 
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SUBROUTINE  QFQUT  ( SOS  -  3.  P.  PFIT.  ERR.  N.  AH) 


SVSCOCOHEYS.  F 
nolist 

SYSC0M7ERRD.  F 
NOLIST 


MNEMONIC  KEY8  FOR  FILE  SYSTEM  (FTN) 
MNEMONIC  COOES  FOR  FILE  SYSTEM  (FTN) 


DIMENSION  SOB  ( t  ) .  S(  1  ),  P  (  1  ) .  PFIT  <  1  ) ,  ERR  {  1  ),  AK(3) 
DIMENSION  NAMEF (6) 

READ  OUTPUT  FILE  NAME, 

CONTINUE 
WRITE!  1.  1001  ) 

READU.  1002 )  NAMEF 

SEARCH  TO  SEE  IF  FILE  ALREADY  EXISTS. 

CALL  SRCH*»(H»nxST.  NAMEF.  12.  1.  u,  ICODEi 
IF ( I CODE  .  EO.  E9FNTF ) 00  TO  30 
WRITE!  1.  1006  ) 

READ!  1.  •)  I  OH 
IF  (  I  OK  EO.  1)00  TO  20 
00  TO  10 
CONTINUE 

OPEN  FILE  FOR  WRITINO. 

call  3RCH*»(K9WRIT,  NAMEF.  12.  1.0,  ICODE) 

IF  (  ICODE  NE  0  >  CALL  ERRPR»(K»NRTN,  ICODE,  O.  0,  0,  0) 

WRITE  RESULTS. 

WRITEO,  1003) 

WRITE (3.  1004) 

WRITEO.  1003) 

DO  30  I-l.N 

URITEO,  1010  )  SOS  (I).S(I!,P(1).PFIT(I).  ERR  (  I ) 
CONTINUE 
WRITEO.  1007) 

WRITEO,  1008) 

WRITEO.  1009)AH(1).  AK<2),  AH  (3! 

CLOSE  OUTPUT  FILE. 

CALL  SRCH49 < H9CL0S.  NAMEF,  13,  1,  0,  ICODE) 

IF ( ICODE  .  NE.  O ) CALL  ERRPR*(H*NRTN,  ICODE.  0- O.  0,  0) 


FORMAT! 'ENTER  OUTPUT  FILE  NAME.  FORMAT  6A2 ' > 

FORMAT <6AP> 

FORMAT!///, '  RESULTS  OF  UFIT') 

FORMAT!//,  '  SNDB  '  SN  '  P  '  f 

*  '  ERR  '  ) 

FORMAT! / ) 

FORMAT! 'OFOUT:  *  EXISTING  FILE  •  IF  OK  TO  MOD  ENTER  1. 
FORMAT!///,  '  FIT  PARAMETERS  ') 

FORMAT!//,  '  P"  1  -  (  1  +  ( A*SN  >  »*B  >  •*<  — C  )  '.  /  / , 

*'  A  3  '  C 

FORMAT (  SOX,  E13  3)  ) 

FORMATO!  IX.  FB.  4,  IX)  ) 


ELSE  0  '  ) 
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SUBROUTINE  QFDFIT  (PD,  SDB,N,P) 

DIMENSION  SC  1 3 > ,  Y(2),  2(2),  W(15),PD(1),P(3), SDD(15) 
DIMENSION  PDD ( 200 ) 

REAL  MV,  MZ 
DO  100  1*1 , N 

S( 1 )=<SDB(  I  )  ) 

PDD  <  I )  =PD  (  I  )  /100.  0 
100  CONTINUE 
C 

C  COMPUTE  SLOPE  BETWEEN  FIRST  TWO  POINTS 

C 

DO  150  1*1, 2 

Y ( I ) “ALOG 1 0 ( -ALOG 1 0 ( 1 -PDD ( I ) ) ) 

150  CONTINUE 

MY* ( Y ( 2 ) - Y ( 1 ) ) / ( S ( 2  > -S ( 1 ) ) 

C 

C  COMPUTE  SLOPE  BETWEEN  LAST  TWO  POINTS 

C 

NM1=N— 1 

DO  200  I =NM 1 ,  N 

Z (  I  >  “-ALOG 10 ( 1 -PDD ( I >  ) 

200  CONTINUE 

MZ* ( Z  < N )— Z (NM1 ) ) / ( S ( N ) — S( NM1 ) ) 

C 

C  COMPUTE  INITIAL  ESTIMATE  OF  K3 

C 

P ( 3 ) =MZ/MY 
C 

C  DISTORT  THE  PD  VALUES 

C 

DO  225  1  =  1,  N 

PDD  (  I  )  =  l-(  (  1  -PDD  (  I  ))■**(  1/P  (3)  >  > 

225  CONTINUE 
C 

C  COMPUTE  LARGE  PD  SLOPE  RATIO 

C 

SRL1 =AL0G10 ( ( ( l-PDD(N) )**(- 1 /P ( 3 ) ) )- 1 ) 

SRL2=A|_OG10( ( ( 1 -PDD ( NM 1 >>**(- 1 /P < 3 ))) -1 > 
SRL3*(-1/P(3>  ) -h-ALOG  1 0 (  l-PDD(N)  ) 

SRL4=( 1/P (3) ) *ALOG 1 0 ( l-PDD(NMl ) > 

SRL=  ( SRH -SRL2  )  /  (  SRL3+5RL4  > 

C 

C  COMPUTE  THE  SMALL  PD  SLOPE  RATIO 

C 

SRS 1 =ALOG 1 0  < ( ( 1-PDD(2  > ** ( - 1 /P ( 3 ) ) ) -1 ) 

SRS2=AL0G  1 0  (  <  (  1-PDD(  1 .)**<- 1 /P ( 3 )))- 1 ) 

SRS3=AL0G 1 0 ( ( -ALOG ( 10.  0) /P(3> )*(AL0G10< 1-PDDC2) ) ) ) 
SRS4=AL0G10( < -ALOG ( 10.  0)/P(3) )*<AL0G10< 1-PDD( 1 ) ) ) ) 
SRS* ( SRS 1 -SRS2 ) / ( SRS3-SRS4 ) 

P (3  >  =P (3 ) * ( SRL/SRS ) 

C 

C  INITIALIZE  COUNTERS  FOR  LINEAR  REGRESSION 

C 

SUM=0 
SUMX=0.  0 
5UMY=0.  0 
5UMX2*0. 0 
SUMX Y=0. 0 
C 

C  COMPUTE  W(I)  AND  PARAMETERS  FOR  A  STRAIGHT  LINE  FIT 

C 

DO  250  1*1,  N 

W(  I  )=AL0G10(  (  (  1-.  0 1  *PD  <  I  )■)**(- 1  /p  ( 3 )  )  >  —  1 ) 
SUM*SUM+i. 0 
SUMX=SUMX+S<  I  > 

SUMY=SUMY+W( I ) 

9UMX2=SUMX2+S( I )*S( I ) 

SUMXY=SUMXY+S< I >  *W ( I ) 

250  CONTINUE 
C 

D=SUM*SUMX2-SUMX*SUMX 

B= ( SUMX2*SUMY-SUMX#SUMXY ) /D 


(0073)  A=(SUMXY*SUM-SUMX*SUMY5 /D 

(0074)  C 

(0075)  C  COMPUTE  K1  AND  K2 

(0076)  C 

(0077)  P<1 )-10. **(B/(10. #A> ) 

(0078)  P (2) =10. *A 

(0079)  C 

(0080)  WRITE  (1,400)  P(1),P(2),P(3> 

(OOSD  400  FORMAT  (  3HK1  =  ,  F6.  4,  2X,  3HK2-,  F6.  4,  2X,  3HK3-.  F6.  4  ) 
(0082)  RETURN 

(0033)  END 
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SUBROUTINE  VBAR 


<0001 ) 

<0002)  C 
<0003)  C 

(0004)  IMPLICIT  DOUBLE  PRECISION  <A-H, O-Z) 

<000  b)  C 

<0006)  DIMENSION  SNR  <  12 ) ,  PDF' <  12 ) ,  TH < 3 ) .  DEL <3 > ,  DINC  < 3 ) ,  W< 3.  12  > 

<0007)  C 

<0008)  COMMON/CCM/SNR,  PDF,  TH,  DEL,  DINC,  W,  M,  II,  JV 

<0009)  C 
<0010)  C 

<0011)  DO  10  1=1, UU 

<0012)  PDF ( I ) *1 .  0+<  < 1 .  DO/TH < 1 ) >*SNR  < I >  >**TH<2) 

<0013)  IF<PDF < I ) .  LT.  1.  DO) PDF ( I ) =1 .  DO 

<0014)  PDF< I >*i.  0/PDF< I )**TH(3) 

(0015)  PDF ( I )*1 .  -PDF < I ) 

<0016)  10  CONTINUE 

<0017)  C 
(0018) 

(0019) 


RETURN 

END 
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1.0 

INTRODUCTION  AND  SUMMARY 

Here  is  a  "quick  and  dirty"  algorithm  for  finding  approximate 
values  for  the  parameters-  K^,  Kt),  and  in  a  Khoury-function  [1]  of  the 
form 

K,  -K 

PD  »  1  -  [1  +  (KXS)  ]  (1) 

where  ?  =  detection  probability  expressed  as  a  fraction 

LJ 

S  =  ratio  of  signal  power  to  background  power. 

Given  a  vector  of  pairs  (PD>  S) ,  the  algorithm  computes  K3  using  the  two 
largest  and  the  two  smallest  values  of  P^.  The  parameters  K^,  and 
consistent  with  this  are  then  found  from  the  full  set  of  data  pairs. 

The  algorithm  might  be  useful  for  computing  starting  values  in 
program  PDFIT  [1]. 

The  parameter  estimates  also  are  useful  "as  is"  for  approximate 
calculations.  For  the  very  few  cases  tried,  the  values  of  ,  K, ,  and  K.^ 
obtained  resulted  in  less  than  about  0.2  dB  error  in  the  required  signal 
to  noise  ratio  obtained  from  the  inverse  of  equation  (1). 


.0 


ALGORITHM  DERIVATION 


The  Khoury-f unction  of  equation  (1)  is  easily  transformed  into 
an  equivalent  logarithmic  Khoury-function  (also  known  as  a  Bucknam- function) 
given  by 

-1/K. 

loglu[(l-PD)  -  1)  -  +  K’  SDB  (2) 

where 

-  K2/10 

K1  ”  K2^°810  H 

SDB  -  10  Iog1Q  S 


(3^ 

(S’) 


The  asymptotic  behavior  of  the  left-hand  side  of  equation  (2)  is 

-1/K. 


l°810[d-PD) 


3  - *  ( -  k  1°*na-h) -  pd*  l-° 


(6a) 


^lo'^V1'  V0-0  (6b> 


This  is  illustrated  in  Figure  1.  Thus,  if  we  were  to  plot  the  quantity 
-log^gd-Pp)  vs.  SDB,  the  slope  of  the  curve  would  approach  K^  •  for 
PQ  near  1.0.  Similarly,  if  we  were  to  plot  the  quantity  log^t-log^d-P^)  ] 
vs.  SDB,  the  slope  would  approach  for  PQ  near  0.  The  ratio  of  the  two 
slopes  is  an  initial  estimate  for  K^. 


This  estimate  of  K^  found  t0  underestimate  the  trie  value  for 
K^  greater  than  1.0,  and  to  overestimate  K^  less  than  1.0.  This  is  caused 
by  the  differences  in  slopes  between  the  true  function  and  its  asymptotic 
forms  as  given  in  equation  (6).  A  correction  factor  can  be  applied  to  the 
initial  K^  by  computing  the  ratio  of  slopes  of  the  asymptotes  and  the  true 
function.  The  new  value  is  then  given  by 

SR. 


K 


3 


K, 


where  SR^  is  the  ratio  of  the  slope  of  the  true  function  to  the  slope  of 
the  asymptotic  function  at  large  P  ,  and  SR  is  the  ratio  of  slopes  at 

U  v> 

small  Py 

The  parameters  and  can  then  he  estimated  by  plotting  the 
left  side  of  aquation  (2)  vs.  SDB  using  this  fitting  a  straight 

line  to  the  data.  Parameters  and  K ^  are  then  found  by  inversion  of 
equations  (3)  and  (4). 


TEST  CASES 


The  algorithm  has  been  tried  on  four  test  cases.  The  teat  cases 
are  the  same  four  as  in  ref  [2].  The  leaat-squarea  fit  parameter  values 
and  the  ”quick-and-dirty"  estimates  are  compared  in  Table  1.  The  errors 
incurred  in  predicting  required  signal-to-noise  ratio  from  the  inverse 
function  using  the  quick-and-dirty  parameters  are  summarized  in  Tables  2-5. 
Note  that  although  there  are  significant  differences  between  the  least- 
squares  and  approximate  parameter  values,  the  difference  in  required  S/N 
is  always  less  than  about  0.2  dB  in  the  four  test  cases. 
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F  D  3 
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0.2294  2 

.521  S.309 

FD(PCT) 
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SNFIT 

ERROR ( D&  ) 

5.012 
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“2.361 

"0.03551 

10.75 

'0 . 972 
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20.93 
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Table  4  - 

■  Swerling  0, 

7  hits ,  ?pA 

-  LO-3 
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PD ( PCT ) 

0.7799 
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Table  5  -  Swelling  2,  4  hit;., 


Table  1  -  Least  Squares  and  Approximate  Parameter 
Values  for  Four  Test  Cases 


Swerling 

Case 

#  of  hits 

PFA 

A 

K1 

K1 

A 

K2 

K2 

a 

K3 

K3 

Maximum 
Error (dB) 

1 

4 

io-6 

.3389 

.3198 

2.181 

2.038 

.3863 

.4210 

.21 

1 

11 

h* 

O 

1 

O' 

.7032 

.6624 

2.248 

2.088 

.  3738 

.4101 

.21 

0 

7 

io~3 

.2294 

.2583 

2.521 

2.563 

8.309 

6.533 

.06 

2 

-4 

io-a 

.1806 

.1790 

3.035 

3.010 

.  7799 

.7925 

.02 

*  approximate  parameter  value 


■  least  squares  parameter  value 


4 . 0  ALGORITHM  COOKBOOK  DESCRIPTION 

Step  1.  For  the  two  smallest  values  of  PqCP^  and  P^,  P^  <  P2)i 
plot  yi  vs  SDB^  where 


yi  ^  lo810^_lo810(1  “  Pi^ 

•  * 

SDBi  -  10  log10(Si)  (8) 


Determine  the  slope  o.^  the  line  connecting  these  cwo  points,  i.e., 

y2  ~  yl 

“y  "  SD«,  -  SDB 


(9) 


Step  2. 
plot  z ^  vs.  SDB^ , 


For  the  two  largest  values  of  PD(P^  and  P^»  P3  <  P^) 1 
where 


-logiQ  (1  -  P^.) 


(10) 


Determine  the  slope  of  the  line  connecting  these  two  points,  i.e. 


(ID 


Step  3.  The  initial  estimate  for  is  found  as 
m 

K,  -  ~  (12) 

J  ID 

y 

Step  4.  Compute  the  correction  factor. 


4a. 


Distort  the  P_  values: 

U  1/K 
P±  -  1  -  (1  -  p1) 


i  »  1,2,  3,  4 


4b.  Compute  large  ?  slope  ratio: 


_  -1/R, 

login[(l-Pj  -  1]  -  log,„[(l-P„) 


-1/K3  log10[l-P4]  +  1/K3  log1(J[l-P3] 


Compute  che  small  P^  slope  ratio: 


5  login[(l-P.) 


log  [iilaiSJ 

8io  1  k3 


-  1]  -  log,„((l-P,) 


'1Ogl0(1"?2)]"1°81Ct‘(iK^)1O810(1^1) 


4d.  Compute  the  correction  factor  and  apply  it  to  K3: 


K  *  _  K  (SRL) 
3  “  *3  (SHS) 


Step  5.  For  all  values  of  P^,  plot  vs.  SDB^,  where 

a  ’1/K3 

Wi  "  lo810[<1"?i)  "  11 

Fit  a  straight  line  through  these  points,  computing  the  slope  and  intercept 


slope  -  a 


intercept  ■  b  -  w(SDB  ■  0) 


Step  6.  The  estimates  for  and  are  found  from  the  slope 
and  intercept  values: 

K_  -  10a 

(Jl.) 

Kx  »  lCP10a; 


An  AFI-  program  implementation  is  listed  in  Figure  2. 
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Figure  2.  APL  Program  Documentation 
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